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LETTER OF TRANSMITTAL. 



Department of the Interior, 

TT. S. GEOLoaicAL Survey, 
Washington, B. C, May 31, 1887. 

Sir : I have the honor to trausinit herewith the results of certain 
investigations, which may be broadly designated as relating to the form 
and position of the sea level. These investigations were begun in part 
previous to my connection with the Geological Survey, but they were 
taken up again in 1885, with your approval, at the request of Mr. G. K. 
Gilbert and Prof. T. 0. Chamberlin, fo* solutions of some special prob- 
lems which arose in their geological researches. The work has been 
prosecuted simultaneously with other lines of office and field work. It 
reached its present form substantially, however, more than a year agoj 
and the principal numerical results of the discussion of Professor Cham- 
berlin's problem are incorporated with his paper on The Driftless 
Area, in the Sixth Annual Keport. The purely mathematical features 
of the paper have been published also in the Annals of Mathematics, 
Kos. 5 and 6, vol. 2, and No. 1, vol. 3. I have delayed offering the 
complete manuscript for publication up to this time in order that I 
might give it a careful revision and check all the more important form- 
ulas by independent processes of derivation. 

The questions treated in this paper are for the most part necessarily 
somewhat mathematical. They are, however, fundamental questions in 
geophysics, and although the mathematical form of presentation has 
been followed throughout, an attempt has been made to state the end 
results and formulas in such a way that they may be understood and 
used with safety by those who may not care to follow the details of the 
analysis. For the benefit of such readers a key to the mathematical 
symbols employed is given iUy addition to the list of contents and gen- 
eral index. 

While the analysis of this investigation was designed especially to 
solve the particular problems of Messrs. Gilbert and Chamberlin, it has 
not been confined to those problems, but has been adapted to the entire 
class of problems to which they belong. It is hoped, therefore, that 
the results of the paper will be of interest and value to geodesists and 
mathematicians as well as to geologists. 

Very respectfully, your obedient servant, 

E. S. Woodward. 
Hon. J. W. Powell, 

Director JJ. 8. Geological Survey. 

(97) 13 



Digitized by VjOOQIC 



Digitized by VjOOQIC 



ON THE FORM AND POSITION OF THE SEA LEVEL 



By E. S. Woodward. 



I. INTRODUCTION. 

1. The problem of the form and dimensions of the sea level surface 
of the earth has been one of peculiar difficulty. The combined efforts 
of the ablest mathematicians of the past two centuries, supplemented 
by the most laborious and costly geodetic measurements have yielded 
us the first approximation only to the complete solution. Fortunately 
this first approximation is exceedingly close. It assigns to the sea level 
a form which differs but slightly from that of an oblate spheroid, whose 
major and minor semi-axes are about 20,926,000 and 20,855,000 English 
feet, respectively. • This spheroid, or reference ellipsoid, as it is some- 
times called, has its minor axis coincident with the earth's axis of rota- 
tion and is usually regarded as sensibly fixed in position and dimen- 
sions. With respect to it the actual sea surface or geoid must be 
imagined to lie partly above and partly below by small but unknown 
amounts, the determination of which, if possible, will constitute a 
second approximation to the figure of the earth. For many if not most 
of the applications of science the reference ellipsoid suffices; the first 
approximation is nearly enough correct. But geodesy, on the one hand, 
has attained such a degree of perfection in precise measurement that 
the discrepancies now brought to light in some of its operations must be 
attributed largely if not chiefly to defects in theory. These discrepan- 
cies must be explained before any considerable advance can be expected 
in our knowledge of the figure of the earth along the present lines of in- 
vestigation. Their true explanation is apparently intimately connected 
with the form of the geoid, and it is to the study of the form, therefore, 
rather than to the determination of the dimensions of the geoid that we 
may look for future progress in geodesy. Geology, on the other hand, 
has raised many questions relative not only to the form, position, and 
fixity of the geoid proper, but also with respect to the allied equipoten- 
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16 FORM AND POSITION OF THE SEA LEVEL. [bull. 4& 

tial surfaces of isolated bodies of water at higher or lower levels. It is 
found, for example, in geological investigations^ that the shore lines of 
extinct seas do not always coincide with existing level lines, but often 
cross them at decided angles, or that the water level lines traced on 
islands in such extinct seas differ in elevation from contemporaneous 
lines traced on their distant shores. Aside from the changes which 
may have been due in these cases to subsidence or upheaval, the ques- 
tion may be raised whether such slopes or differences in elevation rela- 
tive to present level lines may not have been caused by adjacent 
attracting masses, which have since disappeared, like the ice mass of 
the glacial epoch, or in a lake basin by the presence of the water itself. 
Correct and oomplete answers to such questions require a knowledge 
of the existing geoid and of the causes whiXh may have produced secu- 
lar variations in its form and position. 

At present it is Ji)y no means clear how any specially extensive addi- 
tions to our information concerning the more minute features of the sea 
surface are to be obtained. It may even be doubted whether we have 
not reached a practical limit in the first approximatioo ; whether in fact 
the distribution of matter within the earth's surface is not so irfegular 
as to preclude gaining anything more than an empirical formula for the 
deviations of the geoid from the ellipsoid of reference. It seems prob- 
able, however, that the forces producing these deviations have their 
seat in a comparatively thin terrestrial crust resting on a fluid or plastic 
substratum (or nucleus), or that such was the antecedent condition of 
the earth, and that our failure to perceive the relations of the crust to 
the substratum is the chief obstacle to improvement. 

In the absence of a complete rational theory the best evidence which 
analysis can bring to bear on questions pertaining to the geoid is largely 
of a negative character. The effects which would result under certain 
conditions can be computed, but it is not always possible to prove that 
those conditions accord with the actual facts. Investigation must 
proceed to some extent upon doubtful postulates, and computations 
must be made from uncertain data. But notwithstanding this limita- 
tion on the calculations we are about to consider, they will generally 
possess a value in excluding or confirming hypotheses, or in furnishing 
limiting values for the effects of observed causes. 

2. A considerable class of problems concerning the sea level is that 
in which the attracting or disturbing mass is symmetrically disposed 
about a radius of the earth's surface, and is situated on or near the 
surface. As examples of this class we may adduce the two following, 
which led to this investigation : 

(a) Given the dimensions of a lake basin having a circular border. 
When the lake was full of water it left a trace of its surface along the 
border and on an island at the center of the basin. After the water 
had disappeared a line of spirit levels was run between the water trace 
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on the island and that on the border ; what difference in altitude should 
have been found f ^ 

(b) Assuming the accumulation of ice in glacial times to have been 
in the shape of a spherical stratum bounded by a circle, or some sort of 
meniscus symmetrical about an axis, and that the earth's crust did not 
yield under the weight of the ice, what were the resulting distortions 
in the sealevelP 

It will be seen that these problems are essentially the same. They are 
substantially identical also with the problem of the effect of continental 
masses on the sea level, since the continents may be represented, ap- 
proximately at least, as spherical strata having circular borders, or as 
masses of meniscoid shape. 

3. l^e foUowing paper is devoted to the investigation and discussion 
of this class of problems. An attempt has been made to develop the 
theory of their solution so far as is necessary to render practicable the 
numerical evaluation of the characteristic effects of the disturbing mass 
in any special case. In Articles II to XI the theory of the effect of a 
mass in the s&ape of a spherical stratum having a circular border and 
uniform thickness is worked out with considerable detail. The only re- 
strictions imposed on this mass are that its density is uniform, and that 
the ratio of its thickness to the earth's radius may be neglected in com- 
parison with unity. Expressions for the potential of the disturbing mass 
at any point of the disturbed surface are derived in terms of a definite 
integral and in terms of spherical harmonics; and the degree of approx- 
imation of these expressions is investigated. Equations to the dis- 
turbed surface are assigned for the case in which the effect of the re- 
arranged free water is considered, as well as for the case in which that 
effect is neglected. The disturbance in the former case is shown to be 
equal to that in the latter, which is expressed in compact integral form, 
plus a rapidly converging series of additive terms. 

^ This problem was proposed by my colleague, Mr. G. K. Gilbert, to the mathematical 
section of the Washington Philosophical Society, February, 1884. In his geological 
investigations within the area of the Quaternary sea known as Lake Bonneville, Mr. 
Gilbert has found traces of the central portions of the ancient lake surface to be more 
than 100 feet higher than the traces of the contemporaneous surface at its margin. 
A complete consideration of the effects of the causes which might contribute to this 
distortion requires, obviously, a numerical evaluatiop of the depression of the level 
surfaces within the area, due to the removal of the water. 

A solution of the problem was given by the writer before the above-named society 
in March, 1884, and a more complete discussion will be found in sections 39-41. 

2 To what extent the form and position of the sea level may be modified by the mere 
attraction of glacial masses is a question which has been much discussed by geologists. 
It was proposed to the writer by Prof. T. C. Chamberlin, geologist in charge of the 
division of glacial geology, U. 8. Geological Survey. The question is considered at 
some length in sections 42-52, and a preview of the work of the more prominent 
mathematicians who have discussed the problem is given in sections 53-62. The 
principal numerical, results of the writer's investigations are given in Professor 
Chamberlin's paper on The Driftless Area, in the Sixth Annual Report of the U. S. 
Geological Survey. 
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18 FORM AND POSITION OF THE SEA LEVEL. [bull. 48. 

In Article XII the investigation is extended so as to assign the effect 
of any mass of uniform density having a symmetrical distribution about 
a radius of the earth's surface. Particular attention is paid to a class 
of masses whose shapes are assigned by a formula which represents 
fairly well the mass features of the problems (a) and (b) above. 

Under the head of applications, Articles XIII to XVI, the character- 
istic properties of the equipotential surfaces in a lake basin are first 
considered. Then the variations in sea level attributable to continental 
glaciers or ice caps are discussed at some length. The angular radial 
extent of the ice mass is, for the most of the discussion, assumed to be 
380, for the reason that this is the extent of a mass of nearly uniform 
thickness, which would produce the maximum upheaval of the water 
along its border. The external shapes of the various masses, their vol- 
umes, and the distortions of the sea surface attributable to them are 
given in detail. The minimum thicknesses of ice masses of varying 
angular radial extent, requisite to produce average slopes of five feet 
per mile within one degree (69 miles) of their borders, and the extent of 
variation in sea level on the hypothesis of an alternation of glaciation 
at the poles of the earth, are also worked out. 

In the historical note of Article XV, the allied investigations of Arch- 
deacon Pratt, Mr. D. D. Heath, and Sir William Thomson, on the prob- 
lem of glacial submergence are reviewed. The special cases they have 
considered are shown to be easily derived from the general formula of 
Article XII. 

Finally, in Article XVI, a brief discussion of the effect of continental 
masses in distorting the sea level is given. It is shown that according 
as the continents are or are not superficial masses unbalanced in their 
attractive effects, the sea surface must be very irregular or deviate only 
by minute quantities from the ellipsoidal form. It is also shown that 
although a continent whose radial element masses are in a condition 
bordering on hydrostatic equilibrium would produce but slight disturb- 
ances in the position of the sea level, it might nevertheless cause a-con- 
siderable slope of the sea surface, or deflection of the plumb line along 
its border. 

A. THEORY. 
IL MATHEMATICAL STATEMENT OF PROBLEM. 

4. The solution of the general problem outlined in the preceding sec- 
tion depends on the principle of hydrostatics that the potential of the 
forces producing a liquid surface in equilibrium has a constant value 
for all points of that surface. In the case of the earth, if the potential 
of all the attractive forces acting on a unit mass at any point of the sea 
surface be denoted by P, the distance of the point from the earth's axis 
of rotation by I, and the velocity of rotation by a?, the form of the sor- 
fiace will be completely defined by the equation 
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P+iPar^=B, constant. (1) 

The exact value of P in this equation is a complicated function of the 
densities of the element particles of the earth and of the co-ordinates of 
those particles and the attracted point. For the present purposes, how- 
ever, it will be suflQcient to consider P due to a centrobaric sphere of 
equal mass and volume with the earth and concentric with the earth's 
center of gravity. Since we shall only consider relative positions of 
any point on the sea surface, the potential due to centrifugal force, which 
is represented by the second term in (1), may be neglected. 

5. If ao and bo denote the equatorial and polar semi-axes, respectively, 
of the earth's ellipsoid, and Tq the radius of the sphere just referred to, 

^0= ^^75o. (2) 

Using Clarke's values^ of ao and ho we have 

ao=20926062 English feet, 
2^0=20855121 English feet, 
ro=20902394 English feet, 
log ro=7.32020. 

The surface of the sphere thus defined may be regarded as the sur- 
face assumed by a thin film of sea water covering a nucleus whose mass, 
plus the mass of the film, equals the earth's mass. We shall call this 
ideal surface the undisturbed surface. With respect to it the real sur- 
face of the earth lies partly without and partly within ; but so far as 
small relative changes in sea level are concerned it is practically imma-- 
terial whether we refer to the actual closely spheroidal surface or to 
the simpler spherical one. 

6. Let 

If = mass of the earth, 
p^=mean density of earth. 
Then, 

M=i7tro'p^, (3) 

and the equation to the undisturbed surface is 

^=^;rroV.=Oi, (4) 

Ci being a constant. 

Suppose, now, a new mass, 7w, of density p (positive or negative) be 
placed in any fixed position relatively to the undisturbed surface. The 
resulting sea surface will then differ from that defined by (4). To de- 
termine this difference let V be the potential of the disturbing mass m 

* Comparisons of Standards of Length, made at the Ordnance Survey Office, Soath- 
ampton, England, by Capt. A.. R. Clarke, R. E. Published by order of tbe secretary 
of state for war, 1866. 
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at any point of tbe distarbed surface, and let v denote the elevation or 
depression of this point relative to the undisturbed surface. The equa- 
tion to the disturbed surface will then be 

-Ji- + F= Og , a constant. (5) 

n+v 

The difference of this and (4) to terms of the first order inclusive in 

V is 

'0 

whence, putting 

^-(.y-yo)^- (6.) 

Since M/r^^g^ the velocity increment at the earth's surface due to 
the earth's attraction, (6) may be written 

9 

Foin the last two equations is the value of Y when t?=0, or the value 
of F along the line of intersection of the disturbed and undisturbed 
surfaces. If we put 

.+.o=yg=-[. (7) 

This equation represents the elevation of the disturbed surface above 
a spherical surface of equal potential, whose value is 

= 02, 



n— -^0 



since the difference between this and (5) gives (7). 

The constant Fo may be determined from the obvious condition that 
the disturbed and undisturbed surfaces must contain equal volumes. 

It is evident that the equations just derived will hold true if the 
mass m be a part of the earth's mass, so long as the ratio m/M. may be 
neglected relatively to unity. Thus, in the problems we shall consider, 
m may represent the mass of a continent, the deficiency in mass of a 
lake or lake basin, or the ice mass of the glacial epoch. 
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ni. EVALUATION OF POTENTIAL Y— DISTUBBING MASS ON SURFACE 
OP EARTH, OF UNIFORM THICKNESS AND DENSITY, AND WITH 
CIRCULAR BORDER. 

7. The next step in the solation requires the determination of the 
potential Fof the attracting mass for any point of the disturbed sur- 
face, whether without or within the circle which we have assumed to 
define the boundary of the mass. Although the nature of the mass may 
be such as to prevent the water from permeating it freely, the surface 
the water would take if not so restricted is an essential part of the 
disturbed surface. 

In order to derive an expression for F, let the rectangular and polar 
co-ordinates of any point of the attracting mass be defined by the usual 
relations, viz : 

x=r cos d cos A, 
y=r cos 6 sin A, 
z=r sin 6^ 

in which 6 and A correspond to polar distance and longitude, respect- 
ively, the position of the origin being arbitrary. With reference *to 
the same origin, let the co ordinates of the attracted point on the sea 
surface be 

Qc'=r' cos 6' cos A', 

y'=zr' cos 6' sin A', 

z'=zr' sin 6'. 

liD denote the distance between the attracting and attracted points 
and 

cos ^=cos 6 cos ^'+sin 6 sin ff cos (A— A'), (8) 

2)2=r2+r'2-2rr' cos ^=(r-r')2+4rr' sin^ % . (9) 

The volume element of the attracting mass is 

dxdydz=r^dr sin OdddX. 

Hence, if p denote the density of the attractmg mass, a general ex- 
pression for the required potential is ^ 



'="/// 



rHr sin ddddX , .-.^^ 



We must now evaluate this integral. Taking the center of the sphere, 
of reference as the origin of co ordinates, let 

r=ro+tf, 
and (11) 
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in which u and v are small quantities relative to n, v being the same 
of course as defined by equation (6). ^Premising what will be proved in 

Article IV, namely, that we may neglect quantities of the order ~, -, 

^0 Tq 



and (^ — ^ ^ equation (9) gives 

V ^0 y > 



f. 



i)=2ro sin I. (12) 



2 



Prom the first of equations (11) 



dr=zdu^ 
and (13) 

to terms of the order - . 

As to the magnitude of the quantities neglected, it may be remarked 
in passing that r© is in round numbers 21,000,000 feet (see section 5), 
while u and v may be restricted to values less than 100,000 feetj so 
that the fractions neglected will not exceed g^. 

Without loss of generality we may assume the line from which 6 and 
6' are reckoned to pass through the attracted point, and the plane from 
which X and V are reckoned to pass through the attracted point and 
the center of the attracting mass. In this case ^'=0 and yl'=0, and (8) 
gives ^=ft 

By means of this relation and the equivalents in (12) and (13) the in- 
tegral in (10) becomes 

F=rop r r Cdu cos I dddX. (14) 

If the uniform thickness of the attracting mass be denoted by A, the 
limits of u in (14) will be and h. Let the limits of ^, which are ob- 
viously functions of A, be denoted by 6i and 02. The limits of X are 
evidently equal in magnitude but of opposite signs. Hence we have 

F=2rop rduJ\o^-~^de rdX^ArohpJYmn J - sin ^^dX. (15) 

8. To complete the evaluation of (15) it will be convenient to change 
variables. Consider the spherical triangles formed by the attracted and 
attracting points, the center of the attracting mass, and the points in 
which the arc ff cuts the circle bounding the mass. Thus in Figs. 1 and 
2, let P be the attracting and A the attracted points, C the center of 
the attracting mass, and BDQ the bounding circle. Then 

6»=:AP and \=J5AC. 
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Draw CE perpeudicular to AB and pat 



AC=a, 




JBO=/S, 


PE^s, 




B£!=8o, 


CE=p, 




AE=q, 


From either figure 






e:t=q+s, Oi^q- 


-*o» 




whence 








C08 


l^-%- 



6>2=g+«o, 



The right-angled spherical triangles of either figure give 



^^« ^ cos a 

cos q= , 

cos p 



cos ^0= ^^^ » sin »=Sin or sin \. 

COSj? 



Fig. 1. 

The first two of (17) give 

o ^^«2 ^ 1 I COS a 

2cos^^=l+ , 

2 cos jp 



Fio. 2. 



2 8ia«*»=l-22i£, 
2 cos^ ' 



whence 



2 cos ^ sin 



5 ain ^0_|"(COS j? + C08 a)(C08 j?— COS yg)]^ 



cosi> 



From the last of (17) 



d\ 



__ cos pdp 

~'(C08^2^--COS*a')4' 



23 



(16) 



(17) 




(18) 



(19) 



!N"ow, the last of equations (17) and the diagrams show that the limits 
of Pj corresponding to the limits of A, are and «r or and /?, according 
as the attracted point is within or without the circle bounding the at- 
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tracting mass. Hence, if we denote the potentials in the two cases by 
Vi and F2, respectively, the equivalents in (16), (16), (18), and (19) give 

T,=ironp{\ ^^P-'^'^ 'dp, (21) 

'^ j V.cosii— COS ay ^' ' 

a=/3. 

9. The integrals in these equations are in general elliptics of the third 
species. They may be evaluated by the usual processes applicable to 
elliptics, by series, or by mechanical quadrature. 

The integral in (20) presents some apparent difficulty, since the ele- 
ment function is infinite at the upper limit, except when a=/3. Again, 
in case ar=:0, this integral assumes the anomalous form 

10 



yo^y^' 



the value of which is 7t sin ^, as may be easily verified by means of (15), 

(16), and (18). These peculiar features may be removed by the follow- 
ing change of variables, which secures the same constant limits for both 
(20) and (21). 
For brevity put 

fVcosp-cos^Y^fp, (22) 

J rVc08J>-C08/3Ydp. ,!i;, 

I Vcosjp— cos aj ^ • 

Kien, observing that 

sin*^— sin^2 





cosjp— cos/?__ 2 2 




cos »— cos a ^-2 ot _. 2i>' 
^ sm^ _— sin^^ 


put in Ix 






sin|=sin^sin;/i. 


and in J^ 





sin:|r=sin^sin ;/2« 
2 Z 

M08) 
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These give 



a 



dp= 



2Bm~C0Byidyi 



and 



2 sin ^ cos y%dyi 
dp=:-^ j^ » 

(1— sin2^sin2;K2)* 



and the limits for both yx and y^ are and ^. Therefore 



25 



^i 



Ii= 



(l-siu«^8iuVi)* 



(24) 



«<^; 



12= 



2 sin'* ,^ cos* ^'z'^/a 



siu2^ 



sin 



^1 1- ^sinVi I fl-sin^Jsin^^) 



(25) 



10. Some special values ot the integrals (22) to (25) and the corre- 
sponding potentials (20) and (21) are worth deriving. These values are : 

(a) For a point of the disturbed surface at the center of the disturb- 
ing mass, a=0, and (24) gives 



P 



Ji=2sin^ f^^dyi=7rsm/^ 
2J0 ^ 2 



and the corresponding value of the potential is 



Vi=^rohp7t sin 



/? 



(26) 



(27) 



(6) For a point at the border of the disturbing mass, «=/?, and hence 
from (24) and (25) 



Ji= J2=2 arc sin Tsin ^ sin ^'"l "=/?, 



(28) 



(109) 



'y=0 



Digitized byLjOOQlC 



26 FORM AND POSITION OF THE SEA LEVEL. [bull. 48. 

a result which is reached more readily from (22) or (23). The corre- 
sponding value of the potential is 

Fi=y2=4roA/)/?. (29) 

(c) For a point of the disturbed surface ISQo from the center of the 
disturbing mass, a^Tt, and (25) gives 



J2 = 2 




sin« I cos^ yzdyz /^ Yl— cos^^^^i^^j 



l-sin2^ sin« n / 1+cos*^ tanVa 

2 ' f/o 2 ^ 



=;r— 2cos ^ arc tan ftanf cos^tan 



L-V-. -'■.)}: 



=;r('l-cos^^=2;rsin2|, (30) 

and the potential is 

F2=8ro^p;rsin2^. (31) 

{d) Suppose ^=7t; in other words, let the attracting mass cover the 
whole sphere. Then (24) gives 

/i=7r, (32) 

and we have the following well-known approximate value for the po- 
tential of a spherical shell for a point on its interior or exterior surface^ 
viz: 

irohpTT. (33) 

This result follows also from (27) or (29) if we make /3=7r. 

IV. DEGREE OF APPROXIMATION OF THE EXPRESSIONS FOB THE 

POTENTIAL V, 

11. In deriving the expressions (20) and (21) for the potential of the 
disturbing mass, it was assumed that a sulficient degree of approxima- 
tion is attained if quantities of the orders - , - , and upwards are 

neglected. The grounds of this assumption need to be examined with 
some care. For this purpose we shall derive the exact expressions, in 
form at least, for the potential of the disturbing mass at its center, at 
its border, and at 180^ from its center. A comparison of these exact 
values with the approximate values given by (27), (29), and (31) will 
show the order of approximation of (20) and (21). 

^ See eqaation 34. 
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We will first write down the expression for the potential of a com- 
plete spherical shell, of uniform thickness and density, for a point within 
its bounding surfaces. This expression will be useful as a check on 
formulas relating to partial shells. 

Let the radius of the interior surface of the complete shell be ro , the 
uniform thickness A, the uniform density p, and the distance of the at- 
tracted point from the center of the shell n+v. 

Then the potential is^ 

F— 4r hn'r(^l-\- ^ _ '^^ -l ^^^ _ ^^ ^ ^4^ 
'^\^^2{ro+v) 2{ro+v)h^2ro(ro+v) iyro{u+v)hJ ' ^^ 

12. Let the notation be the same. as that used heretofore; i. e., let r 
be the radius- vector of any point of the attracting mass, r' the radius- 
vector of any point of the disturbed surface, p the density, h the thick- 
ness, and /3 the angular extent of the mass,* and 6 the angular distance 
between the attracted and attracting points. For points of the dis- 
turbed surface lying above the undisturbed surface, r will be less than 
r'=ro+i? over the range r© to r', and r will be greater than r' over the 
range r' to Vo+h. Bearing these facts in mind it follows that the exact 
Value of the potential of the disturbing mass for the point where its 
axis pierces the disturbed surface is 



T7- o I 2 7 1 sin MS 



I ^1 Sin Ode 

\ ^'^^ 1 Vf2^y./2_2rr' cos 6 



2p7C 



i. 



=2p7t \ (^4rr'sin2|+(r'-rr-(r'-r)j|;+^ 



/lr«+A 

2p7t I (^^/4rr'sin2.^+(r-r02— (r-r')^-^- 



See Price's Calculus, vol. 3, p. 299. 
(HI) 
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Since 

\ 2 / 

we find by expansion, integration, and redaction, to terms of the first 
order inclusive, 

3/iain § —A 

1+ ^ 



] 4(ro+tJ)sin^ 
F= 4roAp;r sin ^ / ) . (35) 



2 



2(ro+v)sin| 2(ro+t?)A sin 1 1 

If we make /^=;r this expression agrees with (34) to terms of the 
second order. 

Equation (35), it will be observed, differs from (27) by certain terms 
which must be small unless h and v are very large. In one of the most 
important applications discussed in the sequel, ^=10,000 feet, t7=3,000 
feet, and /?=38o. With these values, since Vq is iu round numbers 
21,000,000 feet, the quantity within the parentheses of (35) differs from 
unity by less than yirooT7« I*' >^=60o, which is (see section 23) thb an- 
gular extent of mass required to produce the maximum elevation of the 
disturbed surface at the center of the mass, the quantity within the 
brackets of (35) exceeds unity by less than goVoi using the above values 
of h and V. 

13. Similarly, if the attracted point be at the. border of the attracting 
mass, the exact value of the potential is 

?^'' / J{r-r'f+irr' sm^^ • dX 
(112) 
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in which 

8in2 ^= sin'/gqos^A 
2 1— siu^sinU' 

Introducing this value, the first of the above integrals relative to X 
becomes 




sin2/?cosU+ r^^^ Vcos^y^ 
(r'+r)^^ N l-sinVsin^A ^^• 

Since the numerator of the element function in this integral is greater 
thansin/^cosAandlessthansiny^cosA+^^^cos/?, the value of the in- 
tegral lies between {r'+r)/3 and 

(r'+r)6{ 1+^:111':^^ f' ^^ \ 

V ^'+^ (^ Jo Vi-siu^/^sin^A/ 

Suppose the exact value of this integral is 

(r'+r)/?{l+(?i,^|. 



Likewise, represent the exact value of the second integral relative to 
A..by 

r—r'\ 

r 



(r+r')/j{n-<?.^ 



Then to terms of the first order inclusive the above expression for 
the potential becomes 

The first term of this agrees with (29). 
The quantities Qi and Q2 lie between 




Oand^'^ / ^^ 



a/1— sin2y6fsinU' 



When /?=uj or when the attracting stratum covers a whole heml 
sphere, ^,=§2=0, and (36) becomes 

V=2rohp7r 1 1+ -~-^---J!^ I . (36') 
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This agrees, as it should, to terms of the second order with the half of 
(34). For A=10,000 feet and i?=3,000 feet the quantity in the paren- 
thesis exceeds unity by less than j^o^. When /?<^, ^i and Qawill ex- 

ceed 0, and the degree of approximation of (36) will be somewhat higher 
than that of (36'). 

14. For a point of the- disturbed surface 180° from the center of the 
disturbing mass the exact value of the potential is 

sin ede 



rHr / - , t-^ 



+2rr' cos 6 
(Y^r'-' y/r^'+r'^+2rr' cos P'hj^^- 



When 



(r—r')2<4rr' 008.2 ^ , 
which is the only case we need consider, 

Vr^+r'^+2rr' cos//=2 cos ^ \/rr' ) 1+ -^ ^ 



Srr'cos^ . 




osi '"SI 2'-+---^ 



Srr'cos 

Hence, expanding, integrating, and reducing, there results to terms of 
the first order inclusive, 

F=8.^p.sin^|{l+^^^;, (37) 

the first term of which agrees with (31). Using the values ro=21,000,000 
feet, ^=10,000 feet, and i7=— 1,000 feet (v being here intrinsically nega- 
tive), the factor in the brackets of (37) exceeds unity by 2^0 • 

V. DEVELOPMENT OF POTENTIAL V IN SERIES OF SPHERICAL 

HARMONICS. 

15. The preceding expressions for the potential of the attracting mass, 
namely, F, as defined by equations (20) and (21), are sufficient for most 
of the applications to bo considered in the sequel. They possess the 
obvious advantage of a compact integral form. For some purposes, 
however, it will be desirable to have Y expressed in a series of spheri- 
cal harmonics or Laplace's functions. We may thereby arrive at equa- 
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tions (20) and (21) by a process differing from that followed in Article III 
and establish a useful harmonic development of the elliptic integrals 
Ji and I2 of equations (22) to (25). 

16. Expressions fulfilling the present requirements may be derived 
from equation (10) by expanding D-^ in a series of ascending powers of 

—j and — . Thus from equation (9) we have 

i=^[^o+Pi(p)+P2(jy+ . . .] when r<:r', 

^=-^[Po+P,( J )+P.(0V . . .] when r>r'. 

In these equations Po, Pi, P2, etc., are Laplace's coeflBcients of the 
zero, first, second, etc., order^ respectively. They are functions of the 
angular coordinates ^, ^, A, and V only. 

Taking, now, the center of the sphere of reference as the origin of co- 
ordinates, and supposing the line from which 6 and 6' are reckoned to 
pass through the center of the attracting mass, we shall have for a mass 
of uniform thickness fe, all integrations in (10) independent. For that 
part of the disturbed surface which lies above the undisturbed surface 
r'=ro+i? will fall between the extreme values of r, which are Vo and 
n+A; and hence r<r' for values between ro and n+v^ r>r^ for values 
between n+v and n+h. For that part of the disturbed surface lying 
below the undisturbed surface r>.r'. In both cases the limits of are 
0, and the angular radius /3 of the attracting mass and the limits of X 
are and 2;r. Therefore for that part of the disturbed surface lying 
above the undisturbed surface equation (10) gives 

,a / [P.+P.(r)+P.(^r)^...]«r 

smm / ^^ / [^+^i(7J+^2(0+...]r(Zr5 

Cyo t/r„+w 

Likewise for that part of the disturbed surface which lies below the 
undisturbed surface, equation (10) gives 




[Po+Px(0+P^(0+ . . .ydr; (39) 



r'=ro+i?. 
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For brevity, let the integrals with respect to r in (38) be denoted 3.8 
follows: 

j_{ro±vf—rl 
''"- 3(ro+r) ' 



T _ {.n+vy—ro* 

r_(ro±vf-rl 
''*- 5{v,+vf ' 

Jo= 2 , 

Ji'=(ro+h) {h-v), 

J,'=(r,+t)Moge^, 



(40) 



(41) 



The ambiguous form which J/ assutnes when i=2 receives its proper 
interpretation in the third of (41). 
Similarly, let the integrals with respect to r in (39) be denoted thus: 

do ^ J 

Ji"=(ro+r)/t, 

j,"={n+v)nogJJ±J^, (42) 

^0 



jnJ.U+^ rr„-w-»_(ro+7t)-((-«"| . 

Substituting these equivalents of (40), (41), and (42) in (38) and (39) 
the latter become, respectively, 

V=p(\\uede t\j^-\-J,')Po-\-{Ji+Ji')P,+ . . .]d\, (43) 
Jo Jo 

r=p Tsin Ode C^'[Jo"Po+Ji"Pi+ . . ']dX. (44) 
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Now, it is known from the theory of spherical harmonics^ that the 
general value of P, is 

P,=/;(cos d)fi (cos ^0+ terms multiplied by cosines of multiples of (Jl— A'). 

Since 

r cos i(^-AOe?A=0, 

we have to deal only with the first term of P,. The function of cos & 
or cos d' involved in this first term is defined as follows:^ 

//=:COS 6 J 

The following important relation exists between any three consecu- 
tive values of /, (/i), viz : ^ 

/(;.)==?li:l;//^,(;.)-lzi /_,(;,). (45) 

The remaining integrals in (43) and (44) are therefore of the form 

f,{GOBe)muede I dX==27rfi{M') I Ml^W. 
Jo Jcwfi 

Fi(/S)=f' MM)dM. (46) 



Let 



^ cos/3 

The known value of this integral is* 



^* W=^^[/-i(cos /?)-/,+i(cos/?)] . (47) 

The values of Fin (43) and (44) may now be written thus, replacing 
6^ by a, 6^ or a being the angular distance of the attracted point from 
the center of the attracting mass: 



t=oo 



(J,+J/)Moosa)F,(/3)] (48) 



i=0 

for points of disturbed surface above undisturbed. 

i=00 

r=2pnS^ [J/'MGoa a)F,m (49) 

i=0 

for points of disturbed surface below undisturbed. 

^ See Heine, Handbnch der Kugelfunctionen, Theorie nnd Anwendangen^ Erster 
Band. Second edition, G. Keimer, Berlin, 1878. 

^ Heine, p. 19. ' Heine, p. 91. * Heine, p. 93, 

Bull. 48 3 (117) 
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17. The values of F in the last two equations are exact, but their ap- 
plicability to the problems, we have to consider is limited by the slow 
convergence of the series in the second members. In these problems 
^o» U+% and ro+h are nearly equal, and hence the convergence of the 
series depends almost wholly on the convergence of the functions of a 
and /?. 

If we neglect terms of the order A/ro, v/vq, and upwards, equations 
(48) and (49) become identical and equivalent to (20) and (21), To show 
these facts we expand Ji and J/ of (40) and (41) and obtain 

j^+j/^roh[l~ { {i^l)h-^2iv+Q+l^^ } +0,], (50) 

in which Og represents terms of the second order. Likewise, (42) gives 

J,"=ro/^[l-|(^-l)^+^ 1+O.j . (51) 

Hence, to terms of the first order we have ^ 

and (48) and (49) become 

V=2rohp7t \ V(cosa)^<(/?). (52) 



r=2rohp7rS^ 



Now, if cos 0=cos d cos or+sin ^ sin a cos (A.— A.'), 

CJ'^ ' '^'^' Jo J, ^/1-2C08^+1 

But this integral, as shown by the transformation in section 8, is equiva- 
lent to 

I vcos2>— COS ay ^ I V^osj)— cos/^y 

according as a is less or greater than /?. 

^This inference from (50) and (51) does not appear to be qnite satisfactory. For 
large values of i, J< -f J<' and Ji" are less than r^ ; they are each for i=» . The in- 
fluence of the too great factor r^ in the higher terms of (52) is, however, counter- 
acted by the small factors /< (cos a:)Fj (^) in those terms ; and that the order of approx- 
imation secured in (52) is sufficient is evident from the equivalence of (52) with (20) 
and (21), whose order of approximation has been investigated in Article IV. 

(118) 
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This establishes the equivalence of (52) with (20) and (21), and fur- 
nishes the following development in polar harmonics of the elliptic inte- 
grals in (22) to (25) : 






(53) 



^-/(S|3Siy^i>=l2/'(--)^'(/'): 



In order to show the form of the second, numbers of equations (52) 
and (53) we give below the first four values of /(cos a) and !'<(/?), re- 
spectively. The series of values may be easily extended by means of 
the relations in (45) and (47). 

/o(co8 «)=!, J'o(/?)=l-cosyJ=2 sin2 ^ , 

/i(cos a)=cos a^ ^ J\(/?)=i(l-cos2 /?)= J sin^ /?, 

/^(cos «)= J(3 cos^ flf-1 ), J!j(/?) = J(cos yff-cos^/?) = J sin^/? cos /?, 

/3(cos a:)= J(5 cos^ or— 3 cos a)] i^3(/5)=i(6 cos^ /?— 5 cos* /?— 1).^ 

By means of these values equation (52) may be written thus : 
/-f 2 sin^^+lcos a sin^ /? 

T^2r^pn +i(3 ^^«' ^"^^ «^^' ^ ^^« ^ S . (54) 

]+^(5 cos3 «— 3 cos a) (6 cos^yJ— 6 cos* yJ— 1)\ 

\+ . . . 

VI. EFFECT OF REARRANGED FREE WATER. 

18. In case the disturbing mass is as large as the supposed ice mass 
of the glacial epoch, the attraction of the rearranged free water on 
itself may be appreciable. To determine the exact effect of this attrac- 
tion would be a work of great difficulty even if we had the requisite 
information, namely, an accurate knowledge of the complicated shapes 
of the continents and sea bottom. But we may determine an effect 
which will exceed the probable actual effect by supposing the whole 
surface of the earth covered with a film of water free to assume the 

(119) 
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proper form for equilibrium under the given forces. To fit this ideal 
case formula (6) may be modified in the following manner: 

Let the potential Fin (6) be replaced by F+JF, where JFis the 
potential due to the rearrangement of the water. Likewise replace 
the constant Fo by Fo+ J Fq. Then, if Jt? denote the corresponding 
change in v, equation (6) gives 

i?+Jr=^-A_(F+^F-Fo-JFo). (55) 

Now, v+Av may be expressed by a series of Laplace's functions 
(which are in this case polar harmonics) thus: 

v+Jv=roCo{Zo+Zi+Z2+ . . . ), 

in which Cq is a constant of small numerical value ; and hence, denoting 
the density of sea water by p„ , we have, as shown by Laplace,^ 

jr=.iro%p„7t{Zo+lz,+^Z2 . . . ^^Zi+ . . . ). 

F has already been expressed, equation (54), in a series of Laplace^s 
functions. Denoting these functions in (54) for brevity by Fo, Fi, Fa, 
etc., 

V=2rohp7t(To+Ti+T2+ . . . ). 

Substituting these values of F, JF, and t?+ Jt? in (55) there results, if 
we make the obviously permissible substitution 

the following equation : 

+roCoz/i_3^yh^r. 






)=0. 



+ . . . 

According to the theory of Laplace's functions we must have in this 
equation the sums of the functions of the same order separately equal 

^ M^canique Celeste, Book III, Chap. II, 
(120) 
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to zero, which amounts to placing each line in the equation equal to 
zero. Thus we find 



3 






3,P 



1- 



(2t+l)p. 



whence by summation 




(56) 



This equation expresses the total effect of the disturbing mass in 
altering the sea level, v being the effect which would result if the ocean 
were an infinitely rare fluid, and Jv being the increase over v which 
would result under the assumed conditions. Obviously, v and ^dv may 
be expressed separately. Thus 

Jv=hRf\^—^ -1^\ (57) 

2 PA Aj(2i+l)^-3 §^-3 7 

VII. EVALUATION OF CONSTANTS Fo AND Uq. 

19. We proceed now to determine the constants Vo of equation (6) 
and Uo of equations (56) and (57). 

It has already been stated that these constants are to be determined 
from the condition of equality in volumes contained by the disturbed 
and undisturbed surfaces, a condition whose analytical statement is, to 
terms of the order we neglect, 

2ro';r / v sin adaz=0. 
(121) 
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Sabstituting the value of v from equation (6) in this, there results 

/ Fsinadof— Fo / sin ada=0, 
t/o t/o 

whence 

Fsin ada. (58) 



-\S^ 



(a) The easiest way to evaluate this integral is to substitute for Y its 
value given by equation (54). We get, then, at once 

Fo=4roAp;rsin2 ^, (59) 

since by the theory of spherical harmonics all terms of the series except 
the first vanish in the integration. For the same reason, if we apply 
the condition 



/' 



(t?+ Jt?; sin ada^O 

to equation (56), it will appear that 

l7o=ro=2sin2^. (60) 

(b) The value of Fo may also be found by the following process, which 
is chiefly interesting on account of its complication as compared with 
the process used above. For points within the perimeter of the attract- 
ing mass replace Fin (58) by Fi of (20), and for points outside the perim- 
eter replace F by F2 of (21). Making these substitutions, there results 

Fo=2roftp( / Iimiada+ I Iimiada\. (61) 

Substituting the value of Ji from (24), 




sin' ^ - 8in« - sin* yi 
2 2 ^ 



|sin^;/i5( Sin' 



l-ain«|8in2yi / V 



sin^^sin^^i 



'ty 



2 ^ 2 



1— sin*^sin*vi 
(122) 
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Then the last integral becomes 



— 4cos*^ 






ir 
2=4 COS* ^ 




L<,2_i ti2«i-^2 ^•(l-^)(l+ei)JsinVi 



in which 



^i=sin|, 



and 

sin ^cos;^i 



( 1— sin^^sin* ;Ki ) 

Snbstitnting these limits in the non-logarithmic part of the integral, 
it becomes 

+2 8in» ^ cot* ^ /^L_^ lo" ( l+<^)(l-^i ) 

+ J Sm - cot ^ / ^j^^^ lO^e (i_y(i + j^x • 

Integrating by parts all terms of this expression except the first, we get 



(l-sin»| sinVO* ^, 



y\ 



sin ^ sin ;^i sin ^ 

4 sin* ^( + ^^^ s^^ ( s^*^ ^ '^^^ >^M 

— cot* -^ arc sin Tsin ^ sin ;^i j 



This gives 



/ 



•/7 

Ii sin a(ia=2(sin yJ-yJ cos /J). (62) 

(123) 
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The second integral in (61) becomes by substitution of the value of 
I2 from (25) 

f /i sin ocda 

§ \ cos^ y^dy^ I ^^.^ ±2 

y J (sin^f-«i^'|si'^V»)*(l-8m^§8inV.)* 

a\ \ sin ^ cos y2 sin ^ sin^ ;/2 cos ;/2 V 
=8 sin^l I )cos^ y 2 + .^^ i;,, 

IT 

;/2+^ sin 2;/2-2 arc sin ^sin j »in y2)\ 
r 1— sin^^ sin*;^2 ) sin y2 
=4sin«|/ sin I 

' arc sin (sin ^ sin y2 j 



sm»| 



This gives 

J2 sin ada=:2^7r sin^^- sin /3+/3 cob /i\ . (63) 



I 



The sum of (62) and (63) is 

2;r 8in» | , * 

which, substituted in (61), gives for Fq the same value as (59). 

Vin. EQUATIONS OF DISTURBED SURFACE. 

20. By reference now to equations (3), (6), (20) to (23), and (59) we 
find for the equation of the disturbed surface when the eflfect off the re- 
arranged water is neglected 

.=3;^^(I-siu.f). (64) 

(124) 
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The corresponding expression in polar harmonics is [see equations 
(56), (67), and (60)] 



^'^l^yy^[Yi=A{coBa)F, 



m- (65) 



Under the assumption that the water covers the whole sphere and is 
free to adjust itself as stated in Article VI, the equation to the dis- 
turbed surface is ^ 






=»i(l-«0+«,^ 




(66) 



21. The position of any point of the disturbed surface is thus defined 
by the co-ordinates v and a^ v being the elevation or depression of the 
I>oint relative to the undisturbed spherical surface and a the angular 
distance of the point from the axis of the disturbing mass. Jin (64) 
and (66) is to be computed from (22) or (23) or their equivalents (24) 
and (25), according as the point is within or without the perimeter of 
the disturbing mass. The functions/ (cos a) and J'<(/?) are given by 
(45) and (47), respectively. 

22. The general character of the disturbed surface when the effect of 
the rearranged water is neglected is evident from (64). It is symmetri- 
cal with respect to the axis of the attracting mass. It lies without or 
within the spherical surface of reference according as J is greater or 

less than n sin^ ^. The values of I for the point of the disturbed sur- 

£ace at the center of the attracting mass, for points along the border of 
the mass, and for the point 180^ from the center of the mass, are given 
by equations (26), (28), and (30), respectively. If we denote the corre- 
sponding values of v by the suffixes 1, 2, 3, we get 

.,=3*A(8in|-8in»f), 
a=:0; 

*3=3AA(^2sin»|-sin»£), 
(125) 
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The meaning of these equations may be most readily understood by 
reference to Fig. 3. Thus, if the circle FCH represent (in cross- 




section) the undisturbed sea level surface of the earth, and a stratum of 
matW, as an ice cap, ABDFj be added thereto, the liew sea-level 
surface will assume the form indicated by the dotted line. The values 
of Vi and Vz , as shown in the diagram, are i)ositive, while the value of 
t?3 is negative. If, on the other hand, we suppose the space A'B'iyF 
to be occupied by matter of less density than the average density of 
the earth's crust, as is the case in a lake basin, the disturbed surface 
will fall within the undisturbed surface from F to some line PQ, and 
outside the undisturbed surface from PQ to JT, i. e., t^i and V2 will be 
negative and Vz positive, or, what amounts to the same thing, p in (67) 
will be essentially negative. 

23. It is of interest to inquire what angular extent of mass will pro- 
duce numerical maxima of ^1 , 1^29 and ^3, supposing the thickness h and 
the densities p and p^ constant. By means of the usual criteria it is 
readily found that 



t?]=a maximum for y5=s60o, 

o 

^2=a maximum for sin /?= -, or /J=39o 32', 
t73=a maximum for yS=120o. 



(68) 



24. A glance at equation (66) suffices to show that the effect of the 
free water, if it covers the whole earth, is simply to produce an exag- 
geration of the type of surface defined by (64) and (65). The series in 

(126) 



Digitized by LjOOQIC 



WOODWARD.] EVALtJATlON OF INTEGRALS /i AND U 43 

the third member of (66) expressiug this exaggeration is rapidly con- 
verging on account of the diminishing factor 

3 

2 



Pw 



which is, since -^ is about -^, 



22i+5' 



The essential features of the disturbed surface are, therefore, in any 
case, defined by (64) or its equivalent (65) ; and in most cases the effect 
of the rearranged water may be neglected as unimportant, or as of no 
greater magnitude than the uncertainties inherent in the data for actual 
problems. 

IX. EVALUATION OF TKE DEFINITE INTEGEALS Ii AND Jj. 

25. The equations (67) define the position of the disturbed surface in 
some of its i6ost characteristic points. To define its position at any 
other point we must evaluate the elliptic integral Ii or J2, which per- 
tains to such point. These integrals have already been expressed 
[equation (53)] in a series of polar harmonics, which, if more convergent, 
would sufl&ce for computing Ii or I2. It is easy, however, to derive 
more convergent and convenient series than that of (53), and this is 
the object of the present Article. 

First take Ji of (24). For brevity put 

. a 

sm^ 



and 



sm^ 



0= sm ^ . 



Then by Maclaurin's series, or by the binomial theorem, we readily find 



1=26 f\l-u 



' Ji=26 / (1— Asin^i— BsinVi-CsinVi— • • •)^yi> (69) 

in which 

A= ^O—h"), 



(127) 
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The even powers of sin yi may eacb be expanded in a series of the 
form 

c+d cos 2yi+e cos 4tyi+ . . . 

in which Oy dy e^ etc., are constants. But since 

r 

COS 2nydy=:0j 



£' 



we shall need in these expansions only the values of c. The value of c 
in the expansion of (sin ;/)*" is 



c 



_2n(2n-l)(2n-2) . . . . (n+1) (\^ 



1 . 2 . 3 . . . . w " V2y 



(i 



Applying this formula, and making the integration in (69), there re- 
sults 

ri=6;r(l-.JA-.tJ5-/eO- . . . )• 

Hence if we put^ 

Sr,=i(l-62), 
ff2=^(l-6^)(l+36^), 

^4=tA^(1-6^) (5+9&^+156H356«), 
^5=^5^6(l-«^)(7+12&^+186*+286«+636«), 
^g=^^|3i_^ (X_ j2) (21+3562+506*+706«+10568+2316W). 



(70) 



Ii^hn{X-'g\w^—giw^—g^^— ....), 

. a 

®*^ 2 a 

w= g' 6=sin|? a</3. 

sin § 

This series converges rapidly, except for values of w near unity. In 
a practical application, to be considered presently, wherein /3=3S^ (70) 
gives, using terms up to that in w^\ inclusive, Ji too great by about 5 
per cent, for the case w=zl. But this is the most unfavorable case, 

^ The general yalae of g is 

/+1 .1.3.6... (2n-3) 
[+1. 1.3.5... (2n-5)l . rib* 



^*» - H . 2» . 8« . . . n^ . 23* 



+1.1.3.6 . . . (i!»-9) .1.3. 6 "<»-^><*-^> t«( 

1.2.8 1 

+ • . . 

-1.8.5 (2n-l)6«»» 

(128) * ' * 
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and one moreover for which the exact value of Ji is known from equa- 
tion (28). 

26. By a process entirely similar to that followed above, the expan- 
sion of (25) gives, writing for brevity, 

v=cosec ^, 
fl+J6'(l+v*) 

* +Tei„ft«(35+20»'*+18K^+20K«+3ov») 

+T3f*?5^"'(63+35i'«+30v*+30v8+35»^+63»'W)] 
+ 

If in this expression we put 

^4=T04T+4^»^" +13 1072^ + • • • 
"^S^J St 88 1268 144""!" • • • 
"6^T«2 144+ • • • 

we And 

(71) 

w= ^, 6= sin I , «>/?. 

♦sin I 

This series converges somewhat more rapidly than (70). For the 
case in which /?=38o and for the extreme value w?=l, using terms to 
that in tc?-" inclusive, (71) gives I2 too small by about 3 per cent. 

•A general expression for the nth term within the brackets, beginning with the 
third term, for which w=2, is the following : 

'+1.3 .5. . .(2w— 1) 
+1.3.5. . . (^w— 3) 1 . ny2 

^ ^ \+l . 3 . 5 . . . (2n— 5) 1.3. ^^^T^ v* 
2n(2n-l)(2n-2) . . (n+2) T 1.2 

1-^.2^. 32. . .n^.2-» )+1.3.5. . .(2n-7)l.3.5.<^K!^V6f 



+ 

^+1.3.5. . . (2n— l)y2» 

(129) 
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27. For points near the border of the disturbing mass I^ may be 
expressed by a more rapidly converging series than (71). Thas from 
equation (23) 



Let 



'/( 



^ cos /?— cos a^ 



cos^ 



—cos (Xj 




/5-cos «=2 sin Q^^ sin ('^^=2a. 



2 



cos or) (cos2> — cos a)'^ ( cos j) - cos a)^ 



. . pi>. 



Jo COS©— cos a sin a c;„ Z'^— /^\' 

em (^^- J 

=_sin« /''- ^ -,, 

A (cos2>— cosaf 

/^ (cos J?— COS a)^ / (cos2>--cosa)^ 



dX 
da 






etc.; 
whence 



/ 



(C0S2>— cos a)* siu a da* 



/ (cos ^— COS Of)'* 2 sin^ aV^« eiay' 

etc. 

The integrals in the third and higher terms of the above series are 
thus seen to depend on the integral in the second term. Making the 
requisite differentiations we find, to terms of the third order inclusive, 

T,-ftU ^ . ^^cosa a3(3_2 8in^a) x «^^(^) 

' ^ Vsin a^ 2 sin3 a^ 4 sin^ a "l- '' '/Qg^ ^.^ ("^^^ 



_/ 5a sin /? , 3a* cos a- sin yg , \ 

Vlttsin^a"*" 8 sin* a ■*" * '/ 



smi 

(72) 



a=8inf^sin^. 

J 2 



(130) 
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X. SLOPE OF DISTURBED SURFACE. « 

28. Having derived the requisite formulas for computing the position 
of any point of the disturbed surface, it remains to determine the slope 
of this surface relative to the undisturbed surface. 

Differentiating equation (64) with- respect to a, and dividing, the re- 
sult by the radius of the undisturbed surface n, we get 

dv _ Shp dl ,„ox 

Voda^r^np^ * da* ^ 

This expresses the slope or inclination of the disturbed to the undis- 
turbed surface in a meridian plane through the center of the disturbing 
mass; it also expresses the deflection of the plumb line in the same 
plane. 

In order to apply (73) it is essential to have the general value of 

dl 
d^' 



Since 



■"-(i)' 



^«'"_»«,'> rot '^ 



and bence (70) gives 



ai=-^'^^^J+3g,vA. (74) 



i 



l+4g^w'^y 



Similarly (71) gives 






g=-;rcos^ l+^K.w-'). (75) 



da 2 






29. Equations (74) and (75) will suflBce for the computation otdl/da^ 
except for points near to or at the border of the attracting mass. As 
a approaches equality to /? the above series becoaie less and less con- 
vergent, and finally divergent when £jr=/J or ti7=l. This may be most 

(131) 
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readily seen by differentiating (22) or (23) with respect to a-, and then 
making «=/?. Thus we find 

di 1. . r' dp irw''"~2n' 

Likewise the integrals (24) and (25) become, after differentiating them 
with respect to a and then making a=/?, 

IT 

sec^ yi tan^ yidyi 
(^sec^l+tanVi)*' ^^^ 





sec* yidy^ 
("sec* ^+ta.n* yt Y 



(B) 



Vc*|8ecV*<?ra /"' m^y^U^y^y, 



(^sec» l+tanVa)' / (^sec« | +tan» ;^,y 



This shows the equality of (A) and (B) since (B) is plainly infinite, 
its value being 

[tan yi+faeco i+t&no r^)*"!^ 

30. This failure of equations (74) and (75) for points at the border of 
the attracting mass arises from the fact that the expressions (20) and 
(21), though very approximate for the magnitude of the potential F, 
are not sufficiently general to give an accurate value of dV/da, or the 
attraction in the direction of the arc a for those points. To determine 
the slope of the disturbed surface at the immediate border of the dis- 
turbing mass a special investigation is requisite. 

Since by equations (3) and (6) the slope is expressed by 

dv __3 1_ dV^ /ygv 

r^da 4' VoTtp^ ' Voda^ 

(132) 
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we may derive an expression for the attraction d Y/uia directly. The 
exact expression for tlie horizontal attraction towards the ^xis of the 
mass of any element mass is, nsing the same notation as in Article III, 

4tT^Ar siu« ^ cos* ^ dB cos M\ 
(r-.r')2+4rr'sin*|T 

and the integral of this is dV/r^ia. 
Now, as heretofore, let 



In addition pnt 
so that 



Also let 



whence 



dT]=^dr^ 
77=— t? for r=ro, 
t}:=h—v for r^zTo+h, 



^=2roSin|, 



dS=:ro cos - d0y 



co8|=v^-(i-y- 

Making these sabstitntions and neglecting terms of the order 

i5,^,andr^y, 

the above expression becomes 

^H^dSdrj cos MX 

P a • 

Integrating with respect to 7, and substitating the limits given above, 
there results 

Bull. 48 4 (133) 
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If, now, we suppose the attracted point on the border of the attract- 
ing mass, the limits of S will be 0, and, with sufficient approximation, 
2ro sin /5 cos X=c cos \, say. Integrating with respect to 5^, and substi- 
tuting these limits, we get 

/I \ 1^11 /^ /- , c^cos^A. , ccosA.\ 

1 ^1 1 /^ / - , c* cos^ X , COS A.\ 
+PVCOS XdXloge \yl 1+ ^2 +— ^r~J' 

It remains to integrate these last expressions with respect to A, between 
the limits and ^. An application of the formula for integration by 

parts will readily transform them to elliptics, but since their element 
functions decrease very rapidly from the lower to the upper limit, the 
following process^ will suffice. Consider the integral 

Jr^ iji 1 /^ /T7~?cos^ , ccosA.\ 

in which X is such that f — — y ) °^^y ^ neglected in comparison with 

unity. In the cases we have to consider ( — ^-r i will not exceed yi^ 

\c cos Ay 

if cos A=y^ or A=89o 25' about. Then, since 
, * r I ^ , i? cos^ A. . cos A"! 

the above integral becomes 
/ cos A^A loge -^-^^^ = loge — I COS AdA + I COS ArfA loge coa A 

=Moge — — 1 ) sin A + loge (1+ sin A)4-(sin A— 1) loge cos A. 
But since sin A is very nearly unity, the last expression reduces to 

lo«.(^)-i. 

The error of this integral arising from the use of A instead of ^ as the 
upper limit is less than 

{^ 7t ^\ ^ , f /. . c^cos^A , c cos A \ 
(^_.;iJcosAloge|^l+— ^— -f^j— 1^ 

which, if cos A=ji^ and c cos A/t?=lO, amounts to about y^. 

1 Given ia a somewhat different form by Helmert in Theorieen der hoberen Geo- 
dasie, Vol. II, p. 322. 
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For the entire attractioD, therefore, of the mass for a point on its 
border we get 

i^=-2p{ (.-.)(log,^-l) +.(,og,$-l)} 

Finally, restoring in this last expression the value of o, viz, c= 
2ro sin /?, (76) becomes 

dv 3 ho ( . 8ro sin /? . t? , A— d ^ 1 ,-., 

XI. DISTURBED CENTER OP GRAVITY OP EARTH. 

31. Thns far the disturbed surface has been referred to a spherical 
surface concentric with the earth's center of gravity before the dis- 
turbance arose. In determining the effects of the ice mass in glacial 
times this is the proper surface of reference, since we wish to know the 
distortion of the sea level in those times relative to the sea level in 
preceding and following epochs. If, however, it is desired to consider 
the joint effect in distorting the sea level of existing masses, like the 
continents, on the hypothesis that such masses rest on the surface of a 
centrobaric sphere, a better surface of reference will obviously be the 
disturbed or existing center of gravity of the earth. The use of the 
latter center will require a slight modification of the preceding formulas 
defining the disturbed sea surface. 

To determine the radial displacement of the earth's center of gravity 
due to the addition of such a superficial mass as we have considered, 
it is only necessary to equate the statical moment of that mass to the 
statical moment of the earth's mass, the moment plane being perpen- 
dicular to the axis of the disturbing m^ss at the undisturbed center of 
gravity of the earth. The moment of an elementary ring of angular 
radius /?', measured from the axis of the disturbing mass, is to our 
order of approximation 

2ro^hp7r sin /3' cos /3'd/3' . 
Hence, if a denote the displacement sought and M the earth's mass, 

Ma^n^hpn f 2 sin /3' cos /3'd/3' 



^r^hprt sin^ /?. 

Therefore, by substitution of the value of M given in equation (3), we 
find 

/>« 
(135) 



a^%h^ sin» /5. (78) 
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l^owif the elevation of any point of tUe disturbed suifaee relative to 
the sphere in the new position will be less than its elevation relative to 
the sphere in the former position by an amount whose value to the 
proper degree of approximation is 

or cos ay 

a being, as heretofore, the angular distance of the point from the axis 
of the disturbing mass. That is, if v^ denote what v becomes by the 
change in position at the sphere of reference, 

^' = 1' — (T cosa. 

Hence, by virtue of (64) and (78) we find for the equation of the dis- 
turbed surface when the sphere of reference is concentric with the dis- 
turbed center of gravity of the earth, 

t?'=3^Arl-.sin2 ^ A +COS a cos^ ^\\ . (79) 

And the inclination of the disturbed surface to the surface of refer- 
ence is 

i^ =3A_e_(^5^+^ sin a sinv\ (80) 

XIL EQUATIONS OF DISTURBED SURFACE V^THEN THE DISTURBING 
MASS IS OF VARIABLE THICKNESS. 

32. Throughout the preceding investigations the thickness of the 
attracting mass has been considered uniform. On this account the 
range of application of the formulas derived is somewhat narrow. It 
may be remarked, however, before proceeding to extend the investiga- 
tion to more complex masses, that inasmuch as the data for actual 
problems will be in general more or less uncertain, or to a large extent 
ideal, formulas of a more comprehensive and hence more complex char- 
acter are not specially desirable. Approximate calculations of a rather 
rough sort in some cases will be as good as the data for those calcula- 
tions. The effects assigned by the foregoing equations will be for the 
most part in excess of the probable actual effects, and in so far as com- 
putation can contribute arguments pertinent to observed facts the max- 
imum effects will be most essential. On the other hand it will be 
desirable in some cases to get an idea of the inferior limiting effects. 
The most important of these cases relates to the extent of submergence 
attributable to the ice cap of the glacial epoch. There would seem to 
be little probability of uniform thickness in such a cap. Apparently 
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some sort of regular decrease in thickness (with here and there consid- 
erable though comparatively unimportant deviations) from the center 
to the perimeter of the mass is more probable. Such a law of decrease 
is expressed by the equation 



^'«='=*{'-(E|) 



(81) 



in which h is the thickness along any radial line whose angular distance 
from the axis of the mass is /?, ho is the thickness along the axis, and /3o 
is the angular radius of the perimeter of the mass. In brief, h is a func- 
tion of /?, as stated by the first member of the equation. The exponent 
n must be a positive number and may be for our purposes restricted to 
integer values. In order to determine the effects of masses conform- 
ing to the above, and, in general, any law requiring symmetry of mass 
with respect to a radial axis, we shall devote the present Article to the 
necessary extension of the formulas already derived. 
33. The differential of equation (64) with respect to /J gives 






(82) 



This expresses the elevation of the disturbed surface due to an annu. 
lus of angular radius /?, of angular width d^j and height fe, the density p 
being uniform. If in this equation we make h a function of /?, or write 
h=(p{/3)j and integrate between the proper limits, the result will be the 
elevation of the disturbed surface due to a mass whose thickness con- 
forms to the law expressed by q){/3). Calling, for the sake of distinc- 
tion, v^' the new value of the elevation of the disturbed surface, and 
the proper limits of /?, /3i, and /?2j the result of this integration is 






cp{^)dft. (83) 

This equation assigns the effect of any homogeneous mass whose 
bounding surface is one of revolution abgut a radial axis, subject to 
the restriction that the maximum thickness of the mass may be neg- 
lected in comparison with the radius of the earth. It is obvious, 
however, that the integral in (83) may be impracticably complex for 
some forms of ^(/?). To avoid undue complexity and at the same time 
attain results suitable for our special purposes we shall here confine 
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attention to that form of (p{fi) expressed by eqaation (81). For this 
function we have, considering the whole mass, 



/Si=0, ^2=/So, 
9>(/?)=0 for /? =/?«, 

)=-nV V — V V. 



d<p{fty^ 

sin- I 
Then, observing that J=0 for tf=0, (83) becomes 

The definite integral in this expression depends on and will in gen- 
eral be no less complex than J, which is defined by (22) to (25). An 
examination of (24) and (25) shows that for. points of the distarbed 
sorface within the perimeter of the disturbing mass 

y"l(8in|)'^dsin|=y'j.(8infj-dsinf 

Ji(^sin|-^'*"d[sm|. (85) 

For points of the distarbed surface without the perimeter of the 
mass it is only necessary to replace I in (84) by Jg of (25), or replace 
the limit a in (8c) by y^o- By means of the series (70) and (71), or the 
harmonic series (53), (84) may be evaluated for any point of the dis- 
turbed surface. 

For two points of the disturbed surface, namely, at the center of the 
mass and 180° from that center, (84) yields to direct integration. The 
process of evaluation is as follows : 

The integral in (84) is a function of a^ j3o, and n. Let it be symbolized 
hjf{a/3on). For the two points noted above let this function be dis- 
tinguished by the suffixes 1 and 2, respectively, so that it becomes 

/i(a, /?o, n) for af=0, 

and 

fiiaj/So^n) for a=:7r. 

Then, since from (26), I=:7r sin ^ for af=0, 

/iKA,n)=-^^;rsin|?. (86) 
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Likewise, since from (30), J=2;r siu^ ^ for a=;r, 



/,(a,/Jo,n)=2;rsin^^«-2;r / ( — ^ )(!sio^|. (87) 



(a,/?o,n)=2;rsiii^f -2;r / I ^-j^ \ d 

=1(^4 -2 cos I? - /3o cosec ^ ^ for w=l, 
=|A cos'^ ^« -l\aii2 ^« for w=2, 



When n is more than a few nnits the integral in (87) may be evalu- 
ated by rapidly converging series. Thus, call the required integral JB, 
and let 



Then 



sin ^ = « sin ^^ . 



A / ^ '^°'-^ 







(88) 



Po \ 

5 +. . .y 



+2 ^2 w+4 ^8 n+6 



34. It will be particularly essential for our purposes to evaluate (84) 
for points outside the border of the disturbing mass. The integral re- 
quired is, if we write for brevity 

6=sin^ and &o=sin^. 



" j I2&" ^^ 



Ifow, 1% from the equation preceding (71),. may be written thus: 

I%^7t(jx¥+J2i'+jzh'+ . . . )» 
(139) 
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in whicb ji, J2i etc., are obvious functions of a or y=cosec a. Therefore 
the above integral becomes 

«c r V-''^*^'^ r^-^"^H (89, 

'^ * 1 . « • J. (.1 

\+ • •. ■ 

=;r;S, say. 

35. Let the value of 8 in the last equation for points at the border of 
the disturbing mass where a=/?o be denoted by ^i . Also denote by 
^i"i VS and i?3" the elevations of the disturbed surface at the center of 
the mass, at its border, and 180o from its center. Then equation (84), by 
means of the results in {S6) to (89), gives the following equations analo- 
gous to the group (67) : 

'■"='»^(sTI»'°t-5?5""'f> 
ar=0; 

..-3/^£(«,--^sin^f), (90) 

t?/'=3Ao ^(^2 sin^ 4'-2B ^ sin«^\ 

Pm\ 4 n+2 2/ 

a=;r. 

36. To define the slope of the disturbed surface it is in general nec- 
essary to differentiate (84) with respect to a. The result is of a com- 
l)lex character and subject to discontinuity for points at the border of 
the mass. Foi* practical purposes, however, it will suffice to make use 
of Av"/Aa instead of the differential coefficient, and thus determine 
average slopes over some finite portion of a meridian section of the 
disturbed surface. 

37. In discussing the disturbance of the sea level attributable to the 
ice mass of the glacial epoch, it will be of interest to estimate the effect 
of the rearranged free water. For this purpose we may extend the 
second term of the third member of equation {(oQ) so as to make it as- 
sign the effect of the rearrS^nged water when the mass is of variable 
as well as uniform thickness. The process is strictly analogous to that 
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followed ia deriving (83) and (84) from (64). Thus, confining attention 
to the form of mass defined by (81), and writing, as in section 34, 



-{-(^)] 



and 

we readily fiud from (66) 

9 .vS//<(co8«)«&o- f\(b)b'-'db\ 
' ^-4j\^ (2i+l)|-3 J 

^ few values of Ft(b) derived from (45) and (47) are the following: 

F,{b)=2{y>- b*), 

Fi{b)=2{¥- 36«+ 2b% 
F3{b)=2{¥- 6fr<+10i»- 66»), 
^4(6)=2(6'-106«+306''- 35ft»+ 146"), 
^5(&) =2(6«- 156*+ 706«- 1406«+ 1266'»-42J»). 
From these the corresponding integrals 

/bo ■ ^_j 
F^{b)b db 

can be readily derived. Thas, for example, 

^' b. 



n+4 
70 



F>{b)b''-'db:=2n( ",+® V. (92) 



n+8 
^n+10 • 



42 
"n+12 



h"^ 
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38. A common property of the formulas (84) to (92), both inclusive, 
is worthy of notice. They all refer to ^ mass whose thickness conforms 
to the law (81), namely, 




When n=QO this gives 9>(/?)=ft=fto, or the thickness of the mass is 
uniform. Therefore the formulas (84) to (92) should return to the forms 
applicable to a mass of uniform thickness on making n=x . Such is 
the case. Thus (84) becomes (64), as is readily seen by an application 
of the formula for integration by parts. Likewise, for n=:x equations 
(86), (87), and (89) become (26), (30), and (28), respectively, and the group 
(90) assumes the simpler forms of the group (67). 

B. APPLICATIONS. 

XIII. RELATIVE POSITIONS OF LEVEL OR EQUIPOTENTIAL SURFACES 

IN A LAKE BASIN. 

39. Consider the question stated in section 2, (a), relative to the level 
surfaces in a lake basin. In this case it is required to determine the dif- 
ference in elevation at the center of the basin of two level surfaces 
which intersect along its perimeter. The first two of equations (90) 
give 

V'-t?/'=3Ao^f iSx--^ sin §\ (93) 

This represents the difference in elevation of a level or liquid surface 
at the center and at the border of the basin, p must be understood as 
the excess or defect in density of the liquid relative to the average 
density of the superficial strata of the earth. Thus, if the liquid in 
question be water, 

/)=— (J/>„— 1)=— 1.8, approximately. 

If we differentiate (93), regarding p as variable, the result is 

J(t;/'-^/0=3Ao^(/Si-^sin|>^. (94) 

This expresses the required separation of the two level surfaces in 
question ; i. e., the separation at the center of the basin of two level sur- 
faces which intersect at the border and which are the free surfaces of 
two liquids whose difference in density is J/>. 

To illustrate more fully the meaning of (93) and (94) let ABDE in 
Fig. 4 represent a cross section through the axis of the basin. ACB is 
a circular arc parallel to the section of the sphere of reference. If the 
basin be filled with water the section of the water surface will lie be. 
low ACB as ACB. If the water be removed the corresponding sec- 
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tion of the level surface througb A and B will fall below A CB as AC"B. 
Hence, in the diagram 




Fig. 4. 



40. To get some numerical values, let us assume the following data : 

;to=l,000 feet, /?o=arc of 1°, p„=5-5, 
p = —1.8, ^/>=1? 

this latter being approximate difference in density of air and water. 
Then 

3Ao^=982feet, 
3fto^=545feet. 



The assumed value of /?o is equivalent to about 69 miles measured 
along the earth's surface. 

Now, by means of equations . (86) to (89), we find the following table 
of results corresponding to several values of w, which defines the shape 
of the basin — see equation (81). The results in the fifth column of the 
table express the difference in elevation of a water surface at the center 
and at the border of the basin ; and those in the sixth column express 
the depression of the level surface at the center of the basin consequent 
upon substituting air for water as the attracting mass. 

TahU of valttes showing relative positions of level surfa^s in a lake hasin 140 miles in 
diameter and of 1,000 feel maximum {axial) depth. 



n 


n+1 2 


8i 


n+1 2 


Vi"-V^" 


A(V-»i") 


1 


0. 00436 


0. 00161 


0. 00275 


Feet. 
2.70 


Feet. 
L50 


2 


582 


245 


337 


3.31 


1.84 


3 


654 


298 


356 


3.50 


1.94 


4 


698 


333 


365 


3.58 


1.99 


5 


727 


359 


368 


3.61 


?.01 


6 


748 


379 


369 


3.62 


2.01 


7 


764 


395 


369 


8.62 


2.01 


8 


776 


407 


369 


8.62 


2.01 


9 


786 


418 


368 


3.61 


2.01 


10 


793 


427 


366 


3.59 


1.99 


00 


0.00873 


0.00556 


0. 00317 


3.11 


1.73 
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It will be observed that the differences in the last two columns of the 
above table rise to a maximum between the arguments «=5 and n=9. 
That they should do so is evident from an inspection of equations (89), 
(93), and (94). 

The deflection of the plumb line toward the land along the border of 
fuch a lake, supposing it of uniform depth or n=Qo, would be by 
equation (77) IV or 16", according as the basin is filled with water 
or air, and the slope of the water surface at the immediate border 
would be 0.28 feet per mile. 

41. An interesting inference, which might be drawn from the solution 
of the above problem, is that in triangulating a large lake we should 
expect to see from shore to shore with somewhat less elevations of the 
points of observation than the usual formulas for intervisibility of 
points on the earth's surface would require. This would be a correct 
inference, however, only in case the defect in potential due to the water 
in the lake basin is not offset by an excess in potential due to some 
local or general distribution of matter within the earth's crust. 

XIV. VARIATIONS OF SEA LEVEL ATTRIBUTABLE TO CONTINENTAL 
GLACIERS OR ICE CAPS. 

42. As a second application of the preceding theory we shall investi- 
gate the attractive effects of the ice mass of the glacial epoch, assum- 
ing that the earth's crust did not yield to the pressure of the ice. This 
is the problem of section 2(6). 

This problem in its physical aspects presents two difficulties, the first 
of which has not been alluded to in the foregoing sections, and the 
second only partially considered. 

The first of these difficulties is to account for the enormous quantity 
of water required to form such an ice mass as is supposed to have cov- 
ered our northern hemisphere during the glacial period. This mass has 
been usually estimated as not less than 5,000 feet thick at its center, 
and to have extended 30^ to 90^ from that center. It is generally as- 
sumed to have diminishe'd in thickness with some approximation to 
regularity from the center to the perimeter. The superior limit for this 
shape of mass would be a sheet of uniform thickness, and the inferior 
limit a meniscus increasing slowly in thickness from its perimeter 
towards its axis. Evidently, if the water forming such a mass were 
drawn from the ocean, the latter would undergo a considerable diminu- 
tion in elevation, and this diminution might nearly counterbalance the 
attractive effects of the ice in elevating the water along its border. The 
view advanced by Dr. Groll, however, assumes that there is an alterna- 
tion of glaciation at the poles, the epoch of minimum ice cap at the 
one corresponding to the epoch of maximum ice cap at the other. This, 
granting the sufficiency of ice in the two caps, would make the quantity 
of free water in the ocean substantially constant. In the absence of 
definite information on this point, it must be admitted that considerable 
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ancertaiaty may properly be attribatM to oar compated variations of 
sea level, although the equivalent lowering of the sea, if .the water in 
the ice is drawn therefrom, will be determined for each assumed mass. 

43. The second difficulty arises from the fact that such a large mass 
as we seem compelled to assume for the ice cap would produce an ex- 
tensive rearrangement of the sea water; and we ought therefore, in 
computing the potential at any point of the disturbed surface, to take 
account of this rearrangement. We have shown how to do so in an 
ideal case, which presents eflfects for the elevation or depression and 
slope of the disturbed surface greater than the probable actual effects. 
The actual effects, as we shall indicate, probably lie about midway be- 
tween those assigned by the formulas for the ideal case and those as- 
signed by the formulas which neglect the potential due to the rearranged 
free water. But on account of the difficulty in fixing an exact limit for 
the actual effects, our computed results will be subject to a small range 
of uncertainty, which may be regarded, however, as no greater than 
the inherent uncertainty in the more important data of the problem. 

44. Let us now take for the mean density of the earth p«, for the 
density of ice p, for the thickness of the ice along its axis Aoy And for 
the angular extent of the mass /?o, the following values: 

p^=5.5, fto= 10,000 feet, 

p =1,^ /?o=arc of 380. 

This value of /?o corresponds to about 2,600 miles measured along the 
earth's surface 5 it is also very nearly that angular extent of mass^ of 
uniform thickness, which produces the maximum upheaval ^of water 
along its border. (See section 23.)^ 

45. As we shall compute the effects of masses corresponding to sev- 
eral values of the index w, equation (81), it will be of interest to define 
with some precision the shape of the exterior surface of each mass. To 
do this it will suffice to give the slope of the surface of any mass at 

^ This value for the density of ice is aboat 8 per cent, too great ; but by using it 
in the formulas which do not take account of the potential due to the rearranged 
free water, we shall get results differing only slightly from the probable results. 

2 If the water in the ice cap were drawn wholly from the free sea water the angular 
radius of a mass producing maximum upheaval along its border would be much less 
than the value given by the second of equation {1^)^ for in this case the rise in sea 
level due to the attraction of the cap would be offset partly by the fall due to the 
withdrawal of the water. The difference between the rise and fall just mentioned is 
v^" of (90), minus r of (96), or 

The value of (io which will render this difference a maximum is easily found, but 
the result is of little interest since the corresponding mass would produce effects 
much smaller than the possible effects of continental glaciers. 
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several points in a meridian plane. Tlie differential of equation (81) 
gives 

(95) 



Udft 



nhol '''' 2 \„„, /3 



2r„ 



sin 




cot|. 



This expresses the slope or inclination in a meridian plane of the 
bounding surface of the attracting mass to the spherical surface of ref- 
erence. Using the above values of ho and /?o, and for the radius of the 
earth (see section 4) log ro=7.32020, the following table of values has 
been computed. The slopes are expressed in feet per mile. Several of 
the curves whose slopes are given are delineated (with greatly exagge- 
rated radial scale) in Fig. 5. 



Table showing meridian slopes of exterior hounding surfaces of assumed ice masses, 

equations (81) and (95).] 



ISee 



n 








= 




0° 


6o 


10° 


150 


200 


250 


30O 


350 


38° 


1 


Feet. 
3.88 


Feet. 
3.88 


Feet. 
3.86 


Feet. 
3.85 


Feet. 
3.82 


Feet. 
3.78 


Feet. 
3.75 


Feet. 
3.70 


Feet. 
3.67 


2 


0.00 


1.04 


2.07 


3.08 


4.07 


5.04 


5.96 


6.84 


7.34 


3 


.00 


.21 


.83 


1.85 


3.26 


5.02 


7.10 


9.47 


11.00 


4 


.00 


.03 


.29 


.99 


2.32 


4.45 


7.53 


11.66 


14.67 


5 


.00 


.01 


.10 


.49 


1.55 


3.70 


7.48 


13.46 


18.34 


6 


.00 


.00 


.03 


.24 


.99 


2.95 


7.14 


14.92 


22.01 


7 


.00 


.00 


.01 


.11 


.62 


2.29 


6.62 


16. Og 


25.68 


8 


.00 


.00 


.00 


.05 


.47 


1.73 


6.02 


16.97 


29.34 


9 


.00 


.00 


.00 


.02 


. .23 


1.30 


5.38 


17.64 


33.01 


10 


.00 


.00 


.00 


.01 


.13 


.96 


4.75 


18.10 


36.68 


00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


CO 



The numbers in the above table show that for n=l the slope of the 
bounding surface, as defined by equation (81), is steepest at the axis of 
the mass and decreases slowly from the axis towards the border. For 
values of n greater than unity the bounding surfaces slope up with 
decreasing rapidity from the border to the axis of the mass, the amount 
of slope diminishing to zero at the axis in each case. The features 
here enumerated will hold for any extent of mass, i. e., for any value of 
P(y. It will be observed also that the slope for any value of n is di- 
rectly proportional to the axial thickness fto. Hence the slopes corre- 
sponding to any other thickness than that assumed (10,000 feet) may be 
readily computed from the table. 

46. As to the actual form of the bounding surface of the ice mass 
of the glacial epoch we have no precise information. It is generally 
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assumed, however, that the mass was thickest along its axis, and that 
the thickness decreased with some approach to regularity between the 
axis and the border.^ The slope of ascent along the border has been 
estimated as 10 to 35 feet per mile.^ Nordeuskjold observed a rise of 
7,000 feet in 280 miles from the borcfer of the ice plains of Greenland.' 




Fig. 5. 

Scale for section of 8pliere=jrfTnfoiTTFTT> 

Scale for disturbing mass =7T^W=?TiMfiriny » 

h=h \ 1- ®iBl.^| \ , ^o=10,000 feet. 
( sin« ^/3(f 5 ' 

For curve 1, w=l; for curve 2, w=2; for curve 3, n=10; for curve 4, w=qd . 

This corresponds to an average slope of 25 feet per mile. It seems 
most probable that the slope of such a mass would be steepest at its 
border and diminish gradually towards its center. Whatever may have 
been the actual slopes it is thought that the preceding table affords a 
sufficiently comprehensive variety. Our equations, it is true, will as- 
sign the effects of an indefinite variety of other forms, but in the 
absence of more complete actual data the simple forms whose slopes 
have been computed are considered adequate. 

^The directly opposite view is maintained by W J McGee. See his paper on 
Maximum Synchronous Glaciation, Proceedings of the American Association for the 
Advancement of Science, Vol. xxix, 1880. 

* CroU, Climate and Cosmology, p. 244. 

3 I6i(?., p. 245 
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47. In order to form an idea of the amount of water congealed in tbe 
masses whose shapes have just been defined we will compute the thick, 
nesses of spherical shells of radius n (radius of earth's surface), having 
equal volumes with the ice masses, respectively. As the sea covers 
about three-fourths of the earth's* surface the products of these thick- 
nesses by i will represent approximately the necessary lowering of the 
sea level if the water in the ice cap is drawn from the sea. Since the 
thickness of our assumed ice mass along any radial line at an angular 
distance /3 from its axis'is by equation (81) 



/ 8in-|\ 



the volume of this mass will be expressed by the integral 



2ro'hn I ( 1 J I 



sin /?^/?=4 -^, r^Kn sin» ~. 



The thickness r of a spherical shell of radius u and equal volume 
with the above is given to a suflftcient degree of approximation by the 
relation 



whence 



4roV;r=4-^ro2*o^sin*^, 



(96) 



From this equation with our working values, &o= 10,000 feet and /?©= 
380, we find the following values of r and ^r, corresponding to the 
values of n in the preceding table : 

ToMe showing thicknesses r of spherical shells of equal volume with assumed ice masseSf and 
equivalent lowering of sea level f r. 



n 


T 


It 




Feet. 


Feet 


1 


353 


471 


2 


530 


707 


3 


636 


848 


4 


707 


913 


5 


757 


1009 


6 


796 


1060 


7 


824 


1099 


8 


848 


1181 


9 


867 


1156 


10 


883 


1177 


00 


1060 


1413 
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48. Althoagh, as seen from the last table, the amount of water neo- 
essary to form such an ice cap as we are considering would, if drawn 
from the sea, cause a decided lowering thereof, yet the mass of ice 
would be very small compared with the mass of the earth. Thus, for 
example, the mass of a sheet of uniform thickness 10,000 feet and 38^ 
angular radius is only seofra part of the earth's mass. 

49. To determine the position of the disturbed relative to the undis- 
turbed surface it will be sufficient to compute the elevations of the 
water at the center of the attracting mass, along its border, and 180^ 
from the center, by means of formulas (90). As the slope of the dis- 
turbed surface near the border of the mass is of most importance, we 

'shall compute the elevation of the disturbed surface for a circle of 
points 1° distant from the border, and thereby deduce the average 
slope of the disturbed surface within that distance (69 miles) of the 
border. The separate quantities required in this calculation are given 
in the following table for the same values of n as those used in the two 
preceding tables. It will be remembered that w=qo corresponds to an 
attracting mass of uniform thickness. The values of 8 in the last col- 
umn of the table have all been computed from equation (89), e^ept the 
one for M=ao , which has been derived from (72). 

Table of numerieal values of functions defining position of disturbed surface. 



n 


n+1 2 


n+2 2 


Si 


2(ain»^-B) 


8 


1 


0. 16278 


0.03533 


0. 06078 


0. 01796 


0.05885 


2 


.21704 


.05300 


.09249 


. 02699 


.08946 


3 


.24418 


.06360 


.11218 


. 03243 


.10842 


4 


.26045 


.07066 


. 12568 


.03606 


.12139 


5 


.27131 


.07571 


. 13555 


. 03867 


.13085 


6 


.27906 


. 07950 


. 14310 


.04063 


' .13809 


7 


.28487 


. 08244 


.14908 


. 04216 


.14380 


8 


.28939 


.08480 


.15393 


. 01338 


.14843 


9 


.29301 


.08672 


. 15795 


.04438 


.15227 


10 


.29597- 


.08833 


. 16135 


.04522 


.15550 


00 


0.32657 

«- 


0.10600 


0.21111 


0.05448 


0.19846 



For brevity make the following substitutions : 



^.=»-Ti-t- 



n+2 



sm^ 



A 



jr,= ;Si- 



J* ftin^A 



SlU'' 



w+2 



N,= 8i-8. 
Bull. 48 5 (149) 
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Then from the preceding table we get the following table of values : 



n 


iVTi 


N^ 


N. 


m 


1 


+0. 12745 


+0. 02545 


-0. 01737 


+0. 00193 


2 


+ .16404 


+ .03949 


- .02601 


+ .00303 


3 


+ .18058 


+ .04858 


- .03117 


+ .00376 


4 


+ . 18979 


+ .05502 


- .03460 


+ .00429 


5. 


+ .19560 


+ . 05084 


- .03701 


+ .00470 


6 


+ . 19956 


+ .06360 


- .03887 


H- . oo:oi 


7 


+ .20243 


+ .06664 


- .04028 


+ .00528 


8 


+ .20459 


+ .06913 


- .01142 


+ . 00550 


d 


+ .20629 


+ ..07123 


- .04234 


+ .00588 


10 


+ .20764 


4- .07302 


- .04311 


+ .00585 


00 


+0.21957 


+0. 10311 


-0.05152 . 


+0. 01265 



Now, the factor by which we must multiply Ni, Nz, and N2 to get the 
elevation of the disturbed surface above the undisturbed at the center 
of the mass, along its border, and 180^ from its center is 

• 3*0-^=5454.5 feet: 

and one sixty-ninth part of this factor multiplied by JV4 will give the aver- 
age slope per mile of the disturbed surface within 1^ of the border of 
the ice. Hence we get in the table below the results corresponding 
to the several values of n. The plus sign indicates elevation and the 
minus sign depression of the disturbed relative to the undisturbed 
surface. 

Table showing effects in distorting the sea level of ice caps of the same angular radiuSj 38°, 
<indsame axial thickness^ 10,000 feet, hut of varying external slopes, defined hy equation (81). 





Position of disturbed relative to undis- 
turbed surface. 


Average 

slope per 
mile of dis- 
turbed sur- 
face within 

1° of the 
border of the 

ice mass. 


n 


At centerof 
ice mas. 


Along border 
of ice mass. 


180° from cen- 
ter of ice 
mass. 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
00 


Feet. 
H- 6D5 
+ 895 
-f 985 
+1035 
+1067 
+1088 
+1104 
+11W 
+1125 
+1133 

+1198 


Feet. 
+139 
+215 
+265 
+300 
+326 
+347 
+363 
• +377 
+389 
+398 

+573 


Feet. 
— 95 
—142 
-170 
—189 
-202 
—212 
-220 
-228 
—231 
-235 

—281 


Feet. 

0.15 
.24 
.30 
.34 
.'37 
.40 
.42 
.43 
.45 
.46 

1.00 
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50. As already explained, the numbers in the second, third, and fourth 
columns of the last table assign the position, and those in the fifth col- 
umn the meridian slope of the disturbed surface relative to the undis- 
turbed surface, assuming that the ratio of the density of ice to the earth's 
mean density is ^, and neglecting the eflfect of the rearranged free 
water. If ^ were the correct ratio of the densities, and if the sea cov- 
ered the whole surface of the earth, and were free to arrange itself in con- 
formity with tUe attractive forces, formulas (66) and (91) show that the 
numbers in the second, third, fourth, and fifth columns of the table 
should be increased by about 15, 18, 26, and 18 per cent., respectively, of 
their stated amounts. But the ratio -^f is too great by about 8 per cent., 
so that for this reason the above percentages must be reduced to 7, 10, 
18, and 10, respectively. Again, only three-fourths, at most, of the 
earth's surfp,ce is covered with water, so that the reduced percentages 
must be diminished to three fourths their stated amounts. Moreover, 
the sea could not penetrate the ice mass in such a manner as to produce 
the assumed increase in the potential within its border, and hence the 
reduced percentages must be still further diminished. It is estimated 
that the actual effect can not be greater than two-thirds that which 
would follow if the water were unrestricted. Accordingly, § . f = J the 
above reduced percentages, br 3.6, 5, 9, and 5 per cent., respectively, 
would appear to be liberal allowances for the effect of the free water 
in exaggerating the deviation of the disturbed from the undisturbed 
surface at the points designated, over and above the tabular devia- 
tions. In view of the smallness of these possible increments to our 
computed quantities we need give the question of the effect of the re- 
arranged free water no further consideration. 

The results in the second column of the table are the heights to 
which the water would rise at the center of the ice mass if brought 
within it in any manner, as by a canal, and left free to assume equi- 
librium. If the amount of free water were sufficient it would rise or 
fall to the extent indicated in the third and fourth columns, respectively, 
at points along the border and at the antipodes of the center of the 
mass. 

The slopes given in the fifth column will apply to isolated bodies of 
*free water adjacent to the ice mass, and also to the sea surface in the 
same vicinity, whether there be sufficient water to rise to the height 
indicated in the second column or not.^ 

^Receut observations of the beaches of bodies of water contiguous with the ice 
fields of the glacial epoch indicate uniformly that the water surfaces sloped upwards 
toward the ice. In a careful exploration of the beaches of Lake Agassiz, an extinct 
lake which lay in the valley of the Red River of the North during the glacial epoch, 
Professor Upham has found slopes varying from zero to 1.3 feet per mile. (See Bulle- 
tin No. 39 of the U. S. Geological Survey, on The Upper Beaches and Deltas of the 
Glacial Lake Agassiz, by Warren Upham.) Near the south shore of Lake Ontario, in 
New York, Mr. G. K. Gilbert has observed slopes as great as 5 feet per mile. (See 
Science, Vol. I, p. 222. ) 
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The maximum possible slope of tbe disturbed surface would occur at 
tbe immediate border of a mass of uniform thickness. This slope is by- 
equation (77) for ^=10,000 feet and /?=38o, 1.80 feet per mile, and cor- 
responds to a plumb line deflection of 72''. Since it is uot probable 
that the ice cap presented at its border anything like a vertical wall 
10,000 feet high, we infer ♦that a mass whose maximum or axial thick- 
ness is 10,000 feet would be quite inadequate to produce a slope of 1.8 
feet per mile. 

Again, since all the results in the tabl^ are proportional to the axial 
thickness of the ice, to i)roduce as great an average slope as 2 ieet 
per mile within 1° (69 miles) of the border of a mass having the more 
probable slope defined by the index n^6 to w=10, would require an 
axial thickness in round numbers of 60,000 feet, or 9 J miles. The 
slope at the immediate border in the extremely improbable case of a 
uniform thickness of 50,000 feet would be 5x1.8 feet=9 feet per mile. 

51. To produce slopes as great as 4 or 6 feet per mile for any distance 
from its immediate border an ice cap must have great thickness, which 
implies for any large areal extent a heavy draft on the visible supply 
of water. The minimum thicknesses of ice masses of varying radial 
extent which would produce an average slope of 5 feet per mile within 
lo (69 miles) of their borders are shown m the following table. These 
values are computed on the improbable supposition that the masses 
are of uniform thickness. The volume of each mass is indicated by the 
equivalent lowering of the sea level, or |r deduced from equation (96), 
the value of n being infinite. 

Table showing for ice masses of varying radial extent the minimum thicknesses requisite to 
produce an average slope of 5 feet per mile within 1° of the borders. 



Angular radins 
of mass. 


Minimnm thick- 
ness. 


Equivalent 
lowering of sea 
level, indicat- 
ing volume 
of mass. 


. o 


Feet. 


Feet. 


10 


69, 400 


70J 


20 


52, 500 


2,308 


30 


52, 600 


4,699 


38 


50, COO 


7,065 



For masses having moderate surface slopes near their perimeters the 
axial thicknesses must be about twice as great as the minimum values 
given in the table to produce the same average slope of 5 feet per mile 
within lo of the borders. Thus, for the angular extent 38o, and 
for w=6, say, an average slope of 6 feet per mile would require an 
axial thickness of about 125,000 feet, or 24 miles. This corresponds te 
a lowering of the sea of about 2.4 miles, and if the quantity of free 
water were sufficient, to an elevation of the sea along the border of the 
mass of about 4,000 feet. 
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52. Confining further investigation to our working axial thickness, 
10,000 feet, angular radius 38o, and ratio of densities ^f, we may in- 
quire as to the extent of the variation in sea level at any point of the 
earth's surface on the supposition of an alternation of glaciation at the 
poles. For this purpose it is simply necessary to compute the eleva- 
tion of the disturbed surface at the angular distance a of the point in 
question from the pole or axis of the ice cap and for the point 180O— ty 
and take the difference between the results. Since the maximum effects 
are of most interest, we shall compute the variations in sea level on the 
assumption that the ice cap is of uniform thickness. The results, as 
may be inferred from the second table on p. 66, will not differ materi- 
ally, except for points near the border of the ice, from the results which 
would be derived on the assumption of a sloping mass corresponding 
to the index w=6 to w=10. 

To compute the required integrals we may use forrauljis (70) and (71). 
For /?=38o we find 



log ^1=9.3493, 






logA;i=9.7049, 


Jog ^2=8.7422, 






log A'2=8.8079, 


log ^3=8.3452, 






log A-3=8.3849, 


log 5^4=8.0708, 






log A;4=8.1186, 


log ^5=7.8644, 






log A'5= 7.8904, 


log ^6=7.6987. 






logA'6=7.7033. 


the formula 










sin 
sin 


a 

'2 

'1 





we compute log w for a=lOo, 20^, 30°, etc., to 180^, and then the values 
J, or, more conveniently, the values of I/n sin f? required in (64), readily 

follow. The results are given in tabular form below for circles at inter- 
vals of lOo from either pole. 
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Tahle showing disturbance of sea level atiributahle to an ice cap of 38° angular radius and 
10,000 feet uniform thickness^ and variation in sea level attributable to same mass on the 
hypothesis of an alternation of glaciation at the twopoies. 



Angular distance 

from either pole, 

or a. 


Elevation of dis- 
turbed above un- 
disturbed surface, 
orva. 


Variation in sea level 
from epoch of mini- 
mum to epoch of max- 
imum glaciation, or 
Va—Cm—a. 


o 



Feet. 

+ 1198 


Feet. 
1479 


10 


+1169 


1449 


20 


+ 1076 


1352 


30 


+ 893 


1163 


40 


+ 441 


702 


50 


+ 185 


434 


60 


+ 46 


279 


70 


- 43 


170 


80 


- 106 


80 


90 


-152 


(0 


lOi) 


- 186 


80 


110 


— 213 


170 


120 


- 233 


279 


130 


- 249 


434 


140 


- 261 


702 


150 


- 270 


1163 


160 


- 276 


1352 


170 


- 280 


1419 


180 


- 281 


1479 



A graphical represeutatioa of the results in the last two columns of 
the above table is given in Fig. 6, the section of the disturbed and un- 
disturbed surfaces being so developed that the great circle of the latter 
appears as a straight line, AS, The distances between the two curves 
representing the disturbed surfaces measured at right angles to the 
axis AB indicate the variation in sea level between the epochs of max- 
imum and minimum glaciation at either pole. 

The variation in slope of the sea surface at any point during toe in- 
terval between the extremes of .glaciation will equal the sum of the 
slopes at that point for the two epochs. The maximum variation is, 
however, only slightly greater (about 1 per cent.) than the maximum 
slope already computed, and requires, therefore, no further considera- 
tion. That portion of the disturbed surface having the greatest slope 
is made clearly apparent by the curves in Fig. 6. 




Disturbed 



Suffat 



Fig. 6. 

Horizontal scale, TyruioinTiT ; vertical scale, ^z 
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XV. HISTORICAL NOTE. 

53. The effect of the glacial accamulation ia disturbing the sea level 
has been the subject of considerable discussion. The discussion was 
started apparently by Dr. Groll about 1866 5 and from the widely differ- 
ing results quoted by him in his Climate and Time, Chapter XXIV, one 
might infer that the mathematicians who have attacked the problem 
are completely at loggerheads. There has been, indeed, some diversity 
of opinion as to the proper method'of treating the question, and the 
work of one writer at least is quite erroneous; but part of the appar- 
ent discrepancy in the results quoted by Dr. CroU is drfe also to a cad- 
leal difference in the data used as a basis for calculation. 

The mathematicians whose writings on this subject appear to be 
worthy of especial consideration are Archdeacon Pratt, Mr. D. D. 
Heath, and Sir William Thomson. Their investigations may be found 
in the Philosophical Magazine for 1866; Vols. XXXI and XXXII. 
Those of Archdeacon Pratt are reproduced substantially in his Figure 
of the Earth, fourth edition, pages 236 to 238, while those of Sir Will- 
iam Thomson are given also in Croll's Climate and Time, Chapter XXI 11. 

54. Pratt, — In the following notice of Pratt's investigations we shall 
refer to the fourth edition of The Figure of the Earth (published by 
McMillan & Co., London and New York, 1871). as it is the more recent 
and formal treatise. Although in many respects a valuable text book, 
it is marred by some serious errors, one of which we proceed to point out. 

Pratt's investigation is based on the foil owing erroneous proposition 
(see page 212, Figure of the Earth), namely : 
To prove that the effect of a mass at the earth's surface , whether above or below, 

is to make the sea level rise at any place through a space, , where V is the poten- 

9 
tial of the mass for a point on the disturbed sea level, which is in the same vertical 
line with the place. 

By a process of reasoning to which there appears to be no objection 
he arrives at this equation, namely, 

Y 

r+coust=— , (a) 

y 
in which r is the radius-vector of any point on the disturbed sea sur- 
face, 7 the potential at that point of the disturbing mass, and g is the 
well known velocity increment due to the earth's attraction. The radius 
of the undisturbed surface, supposed spherical, being denoted by a, 
Pratt says : 

Let r=a where F=0 or the horizontal attraction of the mass first becomes appre- 
ciable. 

He thus finds in the above equation 

const =— a, 



and hence 

rise in sea 

(155) 



rise in sea level =r—a-:^—. (ft) 
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But this reasoning is strangely faulty. V is not zero for r=a ; in 
fact it is zero only for points infinitely removed from the surface under 
consideration. He should have reasoned thus : Along the line of inter- 
section of the disturbed and undisturbed surfaces r=a. Call the par- 
ticular value of Faloflg this line Vq. Then we have from (a) 



and hence 



const =—a+ ^, 
9 



r-a=Zi=L^^ (c) 

g 



This agrees with our formula (6)', and the constant Fo is to be deter- 
mined from the condition that the disturbed and undisturbed surfaces 
contain equal volumes. Formula (a) gives results which are too great 
by the constant amount — , i. e., this formula measures heights above 

a spherical surface ~ below the undisturbed surface (see section 5). 

Fo, it will be observed, is never of a lower order than F and can not 
therefore be neglected in comparison with F. 

Although Pratt makes use of the correct principle for determining 
the constant Fo in his article 199, he ignores this principle altogether 
in his article 200, referred to above, and again in his article 213. 

55. In his calculation Pratt assumes a sheet of ice 7,000 feet thick at 
the pole to extend over a whole hemisphere, decreasing in thickness, 
however, as it recedes from the pole in the ratio of the square of the 
cosine of the polar distance. He does not consider the effect of the 
rearranged water. His method of determining the potential F is not sat- 
isfactory. It consists (see articles 90, 91, and 92, Figure of the Earth) 
in a species of mechanical quadrature, by which he computes five special 
values of the attraction of a " hemispherico-spheroidal meniscus,^ whose 
thickness varies according to the law stated above. From these five 
values he derives by the method of indeterminate coefficients a general 
formula for the attraction of the meniscus, and from this formula by 
integration he gets a general formula for the potential of the mass. 
The order of approximation of these formulas is not shown and is not 
evident. From a test we shall apply it is inferred that the approxima- 
tion is so rough as to render the formulas worthless. 

We may readily derive the proper expression for the elevation of the 
sea under the conditions assumed by Pratt from our general equation 
(83). In this equation, if the thickness at the pole be denoted by hoj we 
have 

(15G) 
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Therefore (83) becomes 

- /•- /? 



"-'^.i / 'i"^^'"-' / -^""'^» 





The first term in the second member of this equation is — , and the 

y 

second term is — of equation (c) above. Now, for the ratio of the density 
•f 

of ice to the mean density of the earth -^ Pratt uses i , and hence the 

Pm ^ 

constant by which his results should be diminished, if they were what 
they purport to be, is 

^^7000 t.eet=583 feet. 

Y 
56. We may test the correctness of Pratt's formula for computing — 

by deriving the elevation of the disturbed surface at the center of his 
ice sheet. For this point we have by equation (26) 

J=;r sin I , 



and hence 






The integral in ((?), therefore, is 



and (d) becomes 



*"=tr*»-r2^- W 



This assigns the height of the disturbed above the undisturbed surface 
at the center of the ice mass. Now, 

.^=0.18856, 
15 
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but Pratt's formula (see page 237, Figure of the Earth) gives 0.1189, which 
is only about 63 per cent, of the true value.^ By a fortuitous compen- 
sation of errors, however, Pratt's formula gives a nearly correct result 
for the elevation of the sea at the center of the ice cap, the error in the 
potential for this point being about equal to the constant omitted. His 
formula gives for this elevation 0.1189^0=0.1189x7000 feet=832 feet. 
The correct value is (0.18856— i^)/?o= 737 feet, which is about 13 per 
cent, of itself smaller than Pratt's value. 

57. Heath. — The investigations of Heath may be found in numbers 
covin, CCIX, and OCXIII of the Philosophical Magazine. In num- 
ber covin he develops the theory of the effect on the sea level 
of an ice cap of uniform thickness and ^ny angular radial extent, 
and applies this theory to a numerical example. In number CCIX he 
corrects a blunder by which he was led in his previous paper to the 
conclusion that the ice cap would produce a rise of sea level in both 
hemispheres. His last paper in number OCXIII is chiefly interesting 
as a review and criticism of the work of Croll, Archdeacon Pratt, and 
Sir William Thomson. 

The method followed by Heath is that of Laplace's functions (spher- 
ical harmonics). He takes account of the rearranged water on the 
supposition that it covers the whole sphere and is free to permeate the 
ice mass. He considers masses of uniform thickness only, arguing, 
though not very cogently, that this was the probable form of the ice 
cap. His mathematical processes are correct in principle, but his for- 
mulas defining the position of any point of the disturbed surface are 
rather uninviting to the computer on account of the slow convergence 
of the series used. His series is the same as that in the second mem- 
ber of our equation (66), which we have separated into the definite in- 
tegral and rapidly converging series of the third member of (Q(S). 

1 That the first term in the second member of (d)' is — for the point in question is 

easily proved directly. Thus, to terms of the order we neglect, the potential of the 
mass is 





y_ ^ 27rro sin /? . ^o cos^ /3 . 
2r„8iu| 


r„<J/J 




=2rofeopn- /q 


8 
cos 2 cos2 jSdji 






=4roftp;r^;(^ 






Dividing this by g= 


=J^o|Om'r we get 








9 Pm 15 


y^\if^- 

15 Om 
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68. In his numerical example fleath assumes an ice sheet of uniform 
thickness, 10,000 feet, to extend 30^ in every direction from the pole. 
For the ratio of the density of ice or sea water to the mean density of 
the earth he uses i. That is, in the notation of our equation (66) he 
has 

7^=10,000 feet, 
y^=30o, 

Pm Pm ^' 

Heath gives the position of the disturbed surface at the two poles 
and at points 35^ distant from either pole. He uses fifty terms of his 
series. We will verify one of his results, viz, that assigning the eleva- 
tion of the disturbed surface at the center of the ice cap. With the 
above data formula (66) becomes 

t=» 



At the center of the ice cap 



J=;rsin ^, 



fi{c6s a) = \. 
The numerical values of the first ten terms of 






are, in order, 



+0.02300, +0.00058, 

+ .01203, - .01)008, 

+ .00061, - .00039, 

+ .00358, - .00047, 

+ .00172, - .00041. 

The sum of these is 

+0.04817. 

Therefore the above equation becomes 

v+Jv=95d feet+120 feet = 1,079 feet. 

For the same result Mr. Heath gives 1,078 feet. 

69. In his review Heath first devotes some space to a criticism of the 
views of OroU and the work of Thomson. Croll had apparently held 
the notion that th^ change iu sea level due to an ice cap is essentially 
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equal to the change in position of the center of gravity of the earth 
and cap. This notion is indeed correct for the particular ideal case con- 
sidered by Croll and Thomson, but it is in general quite incorrect. 
Heath points out this fact and concludes that the special case discussed 
by Croll and Thomson affords us an inadequate conception of the actual 
problem. 

With reference to the investigation of Archdeacon Pratt, which hi 
been published in number COVllI of the magazine, Heath says: 

I must confess, with some diffidence, that it appears to me radicaUy erroneous. 

He then proceeds to mention, correctly in the main, but with some 
reservation, the defects of Pratt's method. On one point, however, he 
expresses some doubt as to the adequacy of his own process. He is 
not sure that the order of approximation of his expression for the po- 
tential of the ice cap is sufficient, and thinks this may require closer 
investigation. We have cleared up this point in Article IV. 

60. Thomson. — ^The contribution of Sir William Thomson to the dis- 
cussion of this subject is appended in the form of a note to a paper by 
Croll, On the Physical Cause of the Submergence and Emergence of 
the Land During the Glacial Epoch, published in number CCIX of the 
Philosophical Magazine. The note is brief, but it contains a clear state- 
ment of the essential analytical considerations required in the solution 
of the general problem, and of the decided simplification which results 
when the ice sheet has the special form assumed for the purpose of dis- 
cussion by Croll, namely, that of a hemispherical meniscus, whose thick- 
ness (or density) varies everywhere as the sine of the latitude. In 
analogy with the views of Croll, Thomson devises the following i^eal 
conditions, which are interesting as presenting the mechanical features 
of the problem in their simplest form. He says: 

As an assumption leading to a simple calcalation, let ns suppose the solid earth to 
rise out of the water in a vast number of small flat-topped islands, each bounded by 
a perpendicular cliff, and let the proportion of the water area to the whole be equal 
in all parts. Let all of these islands in one hemisphere be covered with ice, of thick- 
ness according to the law assumed by Mr. Croll, that is, varying in simple proportion 
of the sine of the latitude. Let this ice be removed from the first hemisphere and 
similarly distributed over the islands in the second. 

Thomson gives no analysis, but continues: 

By working out according to Mr. CrolVs directions, it is easily found that the 
change in sea level which this will produce will consist in a sinking in the first hemi- 
sphere and rising in the second through heights varying according to the same law 
(that is, simple proportionality to sines of latitudes), and amounting at each pole to 

where * denotes the thickness of the ice crust at the pole, i the ratio of the density 
of ice, and w that of sea water to the earth's mean density, and oo the ratio of the 
area of ocean to the whole surface. 

61. We may readily get Thomson's result from equation (66). Thus, 
the change in position of any point of the sea surface will be, if we 
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represent the thickness of the ice at any point whose polar distance is 
fi by *o cos /?, 



/ 



^fi!^»,c«s/J,J^= 



But from equation (46) 

^W)=y;(cosy^)sin/?5 

and by the theory of Laplace's functions 

ft 
f, (cos 13) cos ^ sin y^(!/?=0, 



/■ 



except when t=l. In this case, since/ (cos /?)=cos /?, the last integral 
becomes 



/ 



cos^/? sin /?dy?=|. 

"^ 



Therefore, observing that/ (cos a')=cos a, <^ being the polar distance 
of any point of the sea surface, we get 



/ 



.A. 



r "_K-3».^^,^=^^ 



COS «. {e) 



Now, since the ratio of the area of the ocean to the whole surface of 
the earth is assumed in Thomson's problem to be o?, we must replace f> 
in (6) by (I— a?)/) and p„ by Qop„. Making these substitutions, and put- 
ting ar=0, the second member of (e) becomes 

^. (/> 

This is Thomson's result. 

Knowing the fact expressed by (e)j namely, that the transfer of such 
a meniscus as we are considering from one hemisphere to the opposite 
one would change the sea level at any point by an amount proportional 
to the cosine of the polar distance of that point, we may get the result 
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(/) ^y equating the sum of the moments of the transferred ice and 
water to the moment of the whole earth. In other words, this partic- 
ular case requires only the amount of shifting of the earth's center of 
gravity. Calling this amount (Ti , and taking a plane perpendicular to 
the axis of the meniscus at the disturbed center of gravity of the earth 
as moment plane, the equation of moments is 



(l-a,)py^ 



27tro sin 6 . Todd . ^o cos 6 . n cos + 

r^" 4 

Gop„ j 2;rro sin 6 . rdO . cTi cos . r cos ^ = 3 TtVo^p^ffi. 
This gives 

0-1 = ^, 

which is the same as (/). 
In his numerical example Thomson takes 

^0=6,000 feet, 

2 
^=3, 

^1, and A, 2^ 
p^ 6' p^ 11 

These data give 379 feet as the change in sea level at the pole during 
the interval between the epochs of minimum and maximum glaciation. 
This is the greatest change in sea level that could occur under the as- 
sumed conditions. 

62, It is to be observed that the numerical results of Pratt(corrected), 
Heath, and Thomson are not directly comparable with each other nor 
with the results we have computed in Article XIV, since they are all 
based on different data. They represent effects due to causes of the 
same kind, but of widely differing magnitudes. No statement of these 
results would be intelligible without an accompanying statement of the 
data on which they rest. To show clearly how widely different writers 
differ in their data the latter have been collected in a tabular form be- 
low. It will be remembered that Pratt took no account of the effect of 
the rearranged water :^ 

^ One of tlie most important contributions to the disonssion of the eflfects of conti- 
nental glaciers on the sea level has appeared since the mannscript of this paper was 
placed in the printer^s hands, viz : Die (Jeoiddeformationen der Eiszeit, von Erich 
von Drygalski, Dr. Phil. W. Pormetter, Berlin, 1887. 
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Table^ showing data used by different authors in discussing theprohlem of glacial submer- 
gence. 



Author. 


Ratio of 

density of 

ice to mean 

density of 

earth. 

P,n 


Ratio of 

density of 

sea -water to 

mean density 

of earth. 

put 

pm 


Angular 
radius ^f 
ice mass. 


Thickness of ice at 

angular distance /3 from 

its axis. 


Pratt 






o 
90 
30 
90 

38 


Feet. 
7,000 cos*^ 
10, 000 uniform. 
6, 000 cos ^. 


Heath 

Thomson 

^Woodward 


h 
ft. 


1 



XVI. VABIATIONS IN SEA LEVEL ATTRIBUTABLE TO CONTINENTAL 

MASSES. 

63. To illustrate the use of the theory developed in determiuing the 
disturbance of the sea level attributable to a continental mass we shall 
consider, in addition to some observations on the general features of 
the subject, the special case presented by the largest of the continents, 
namely, Europe and Asia. 

In contemplating this problem it is important to distinguish two ex- 
treme hypotheses relative to the nature of the earth's crust. On the 
one hand, we may suppose that a continent is simply a superficial ag- 
gregation of matter, which, if removed, would leave a sensibly centro- 
baric spheroid j in other words, the presence of a continent does not 
imply that beneath it the earth's crust is any less dense than beneath 
the ocean. On the other hand, we may suppose that the several radial 
element prisms of the earth's crust are in a state approximating to hy- 
drostatic equilibrium, and hence, that the mere existence of a continent 
implies a defect of density in the strata beneath it. According as we 
proceed from the one hypothesis or the other we shall arrive at widely 
differing results, which may be regarded, however, as the limits between 
which the facts lie. 

64. Assuming, in accordance with the first hypothesis, that the forma, 
tion of the continent in question involved a transfer of the earth's 
center of gravity towards the center of the continent, it will be of in- 
terest to compute the position of the disturbed sea surface with respect 
to spherical surfaces of the same radius and concentric about the origi- 
nal and disturbed centers of gravity, respectively. For this purpose 
we may use formulas (64) and (79). These are based on the assumption 
that the attracting mass is of uniform thickness. A more reasonable 
assumption is, perhaps, that the continents slope up rather rapidly, but 
not abruptly, from the sea shore, attaining their average height at no 
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great distance inland; and more accurate knowledge than we now 
possess might render it desirable to use some form of the more general 
equation (83), which takes account of variations in thickness of the 
attracting mass. As we can not expect, however, at present, to repre- 
sent the actual shape of a continent very closely, and as we shall not 
attempt to estimate the attraction of the rearranged free water on itself, 
it will be best to confine attention to masses of uniform thickness, which 
give effects greater in general than the probable actual effects. 

65. For the angular radial extent and the relation of thickness and 
densities for the continent of Europe and Asia we take the following 
data substantially as they are given by Helmert in his Geodasie^, Part 
II, pp. 313, 314, viz : 

/?=arc of 380, 
37t/>//o„=4,000 meters 
=13,124 feet. 

The last expression is arrived at by taking for the average depth of 
the sea 3,438 meters, and for the average elevation of the continent 440 
meters, the density of the continental mass being assumed to be half 
the earth's mean density, or 2.8, and that of sea water 1. Thus we 
have a mass 3,438 meters thick, whose effective density is 2.8—1=1.8, 
and an additional mass 440 meters thick of density 2.8. These two are 
equivalent, so far as their potential to terms of the order we neglect is 
concerned, to a single mass of density 1.8 (or JpJ and 4,120 meters 
thickness. Hence, in round numbers, the relation above. 

The value for the angular radius /? of the continent is equivalent to 
about 2,600 miles, measured along the surface of the earth. It is the 
same radius assumed in Article XIV for the ice mass. This value is 
also very nearly that angular extent which a continent of uniform 
thickness must have to produce the maximum upheaval of water 
along its border (see section 23). 

The position of the disturbed surface relative to the two spherical 
surfaces concentric with the original and disturbed centers of gravity, 
respectively, will be given at intervals of 10^ from the center of 
the continent to the point opposite, or 180^ from that center. The 

requisite values of I /n sin ^ for this purpose have been computed 

from formulas (70) and (71), using the logarithms of gfi, g^^ etc., /ti, Ih^ 
etc., given on page 69. To determine the slope of the disturbed surface 
with respect to the spherical surfaces of reference (or the deflections of 
the plumb line) at intervals of 10^ from the center of the continent, 
use has been made of equations (73) and (80), the differential coefficients 
dJ/flJa being computed from (74) and (76). Omitting the details of 
the computation, the nature of which is readily apparent from the 

* For lull title see page 86. 
(164) 
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equations referred to, the results are embodied in the following table. 
The first column of the table gives angular distances increasing by in- 
crements of IQo along a great circle from the center of the continent. 
The second and third columns give for each angular distance the elevation 
or depression of the disturbed surface relatively to the spherical sur- 
faces concentric with the undisturbed and disturbed centers of gravity 
of the earth, respectively, elevations being indicated by the plus sigii 
and depressions by the minus sign. The fourth and fifth columns give 
for each angular distance the deflections of the plumb line or inclina- 
tions of the disturbed surface to the spherical surfaces of reference. 
The signs of these deflections are minus or plus according as the angle 
between the plumb line or normal at any point of the disturbed surface 
tod the axis of the continent is greater or less than the angle between 
the radius vector of the same point and the axis of the continent (see 
Fig. 7). 



^Tahle showing the superior limiting effects in disturbing the sea level attributable to the 

continent of Europe and Asia, 



Angular 

distance from 

center of 

continent. 


Elevation or depression of dis- 
turbed surface with respect to- 


Inclination of disturbed sur- 
face with reference to- 


Spherical sur- 
face concentric 

with undis- 
turbed center 

of gravity 
of earth. 


Spherical sur- 
face concentric 

with dis- 
turbed center 
of gravity 

of earth. 


Spherical sur- 
face concentric 

with undis- 
turbed center 

of gravity 
of earth. 


spherical sur- 
face concentric 

with dis- 
turbed center 
of gravity 

of earth. 


o 


Feet 


Feet 


/, 


// 


00 


4-2,881 


+1,637 


- 0.0 


- 0.0 


10 


-i-2,812 


+1, 587 


- 7.7 


- 5.6 


20 


+2,588 


+1,419 


-17.7 


-13.5 


30 


+2,148 


+1,072 


-33.9 


-27.7 


40 


+|;099 


+ 146 


-51.3 


-43.5 


50 


+ 444 


- 356 


-24.2 


-14.8 


60 


+ 111 


- 511 


-14.7 


- 4.1 


70 


^ 104 


— 529 


-10.0 


+ 1.6 


80 


- i235 


- 471 


- 7.2 


+ 4.9 


90 


- 365 


- 365 


- 5.4 


+ 6.9 


100 


- 448 


- 232 


— 4,1 


+ 8.0 


110 


- 512 


- 87 


- 3:2 


+ 8.4 


120^ 


- 661 


+ 61 


- 2.4 


+ 8.2 


130 


- 600 


+ 200 


- 1.9 


+ 7.5 


140 


- 629 


+ 324 


- 1.4 


+ 6.5 


150 


- 650 


+ 426 


- 1.0 


+ 5.1 


160 


- 665 


+ 504 


- 0.7 


+ 8.5 


170 


- 673 


+ 552 


- 0.3 


+ 1.8 


180 


- 676 


+ 568 


- O.i) 


+ 0.0 



Ball. 48 6 
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66. The relations of the disturbed surface to the spherical surfaces 
of reference are shown in Fig. 7. The full-line circles in this diagram 
represent great circles of the spheres of reference on ascale of 2^0,000000 • 
The position of the continent is indicated by the arc AB of the circle, 
whose center G is the undisturbed center of gravity of the earth. The 
disturbed center of gravity, or the center of the other circle, is at C. 
The position of the disturbed surface is indicated by the dotted line. 
The radial distances of this line from the circles of reference .and the 
distance between the centers G and G' are exaggerated 2,500 times. 
QFR represents a normal to the disturbed surface at P. 

67. For the elevation of the disturbed surf^^ce at the immediate border 
of the continent we find from the second of equations (G7) and from 
(79) 1,380 feet or 400 feet, according as the sphere of reference is con- 
centric with the undisturbed or with the disturbed center of gravity of 
the earth. 

68. It should be remarked in this connection that the results given 
in the preceding table, and, indeed, all those in this section pertaining 
to the continent of Europe and Asia, are such as would exist under the 
assumed conditions if there was no other continent. The complete 
problem requires the determination of the resultant action of all the 
continents at any point of the sea surface. For a method of determin- 
ing this resultant, when the components due to the several continents 
are known, the reader may be referred to Helraert's Geodasie. 

(IC6) 
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69. To determine the deflection of the plumb line at the level of the 
sea along the border of the continent, we observe that according to 
our assumed data the portion of the continental mass lying below the 
sea level is about 11,000 feet thick and has an effective density of 1.8, 
while the portion above the sea level is about 1,400 feet thick and has 
an effective density of 2.8. We may compute the deflection of the 
plumb line due to each of these portions by means of equation (77) and 
add the results to get the total deflection. For the lower mass we have 
for use in (77) 

^=11,000 feet, />=1 8, i?=ll,000 feet ; 

and for the upper mass 

/i= 1,400 feet, . /a=2.8, i?=:0 feet. 

Hence the deflections due to the upper and lower masses are 138'' and 
34'', respectively, and their sum is 172". This deflection is relative to 
a radius drawn to the undisturbed center of gravity of the earth. 
The deflection relative to a radius drawn to the disturbed center of 
gravity is by equation (80) 8" less, or 164". 

70. Under the conditions of the second hypothesis, which supposes 
the several radial element prisms of the earth's crust in a state border- 
ing on hydrostatic equilibrium, it is evident that the disturbances of 
the sea level attributable to a continental mass must be of a low order. 
They must, in fact, be confined to terms of no higher order than those 
which have been njeglected in our equations defining the position of the 
disturbed sea surface. The precise evaluation of these terms would offer 
difficulties practically insuperable in all cases, except those which present 
the simplest arrangement of densities in the element prisms of the earth's 
crust. We may form a sufficiently definite idea of their smallness, 
however, by considering the ideal question of the effect on the sea level 
of the radial transfer of a stratum of the earth's crust from some posi- 
tion below to some position above the sea surface. 

Let the stratum considered be a portion of a spherical shell; let its 
border be circular and of angular radius /?. The effect on the sea level 
will obviously be greatest at the axis of the stratum. Hence, we only 
need to derive the change in potential at the sea level and at this axis 
due to the transferred stratum. Suppose the stratum raised to a height 
hi above sea level. Let its uniform thickness in this position be ti and 
its uniform density Pi. Then, r being the radius vector of any element 
of the mass and r' the radius vector of the sea surface, the potential 
of the stratum at its axis and at sea level will be 



2fH7r / f^dr / sm g d^ _ _ 

t/r'+A. J^ \/r2+r'2-2rr'cos^ 

^or'+Ai+e, 

2px7r I r \/rH Y'^^-rr' cos /:^-(r-r') 1 ?^ . 
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This to terms of the first order inclusive is 

If we suppose the stratum in its original position is below the sea 
surface a distance Ihj and that its uniform thickness and density in this 
position are tz and p2, respectively, the resulting potential at the same 
point of the sea surface is 

2 

Since a^i=P2^2; i. e., since the mass of the stratum in the two posi- 
tions is the same, though its thickness and density may vary, the differ- 
ence of (97) and (98) is 

irWi^ sin |/3/i+3f.+6/..+C/.3,^^ .^^^ 



/ 3fi+St2+6hi+(}h2 fi-f24-2/<i-2fe.A 
V ^^' 4r^sinf / 



This expression shows that a stratum of the earth's crust might be 
transferred through a considerable distance radially without materially 
affecting the sea leveh 

Thus, for example, suppose 

/*i=/i2=10,000 feet, f,=f2=5,000 feet. 

Then the fraction in the parenthesis of (99) becomes j^. In other 
words, the potential at the sea surface under the center of the stratum 
would be only ^ as great as that due to an uncompensated external 
stratum or internal vacuity of equal effective mass. 
In the case of the continent of Europe and Asia we have (section 69) 

/«i=0, fi=l,400 feet, • pi=2.S. 

These in (97) give 



4r';rsin ^ 1400x2.8 



('+2oioO 60000 Sin §)■ ^''^ 



Likewise for use in (98) we have 

^2=0, ^2=11,000 feet, pi=1.8. 

These give 

Ar^Tt sin ^ 



|„000..,8(^>-,4,. > y 



(t) 



(ICS) 
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The sum of (a) and (6) reduced to unit density is 
ir'Tt sin d 23720/^1-^ L 

'"^ ftOOft Rill 



)/t 1 . 1 \ 

V ^^ 9000 siu§/ <") 



This is the potential to terms of the first order, inclusive of the whole 
continent at its center and at sea level, on the supposition that there 
is no defect in density of the underlying strata. But if we suppose a 
defect in density of 0.237 uniformly distributed through a depth of 
100.000 feet below the sea level, we must diminish the result (c) by 



The remainder is 



f 23720 A 1 1 Y 

2 .\^ ^»0 840 sin jy 

J y6i)J 900 sin ^^V 



^r'Tt^ 



or about ^ of the potential due to the uncompensated mass. Under 
these conditions of compensation, therefore, the elevation of the sea at 
the center of the continent of Europe and Asia would be only about 10 
feet, whereas the first hypothesis would Require an elevation of about 
2,900 feet (see section 66.) 

71. Notwithstanding the feeble effect a continental mass, whose radial 
elements are in a state approximating to hydrostatic equilibrium, would 
have in elevating or depressing the sea surface, it is conceivable that a 
<K)nsiderable deflection of the plumb line might be produced by such a 
mass along its border. If, for example, we suppose the visible mass of 
the continent of Europe and Asia to be compensated by a defect in 
density of 0.18 uniformly distributed through a depth of 135,000 feet 
below sea level, the deflection of the plumb line along the border of this 
continent would be diminished by 117", leaving still a deflection (section 
69) of 172"— 117''=55" relative to the undisturbed sphere of reference. 
This deflection, 55", it should be observed, is the maximum possible 
value under the assumed conditions. If the degree of compensation as- 
sumed actually exists, it is probable that the real maximum deflection 
of the plumb line is much less than 55", since our calculation premises 
a vertical coast wall for the continent, whereas it presents for the most 
part only a moderately steep slope along the sea shore. 

XVII. LIST OF AUTHORS CONSULTED. 

72. The authors whose works have been specially consulted in the 
preparation of .the preceding pages, and the full titles of their works, 
dates of publication, etc., are named below : 

Brans, Dr. Heinricli. Die Figur der Erde. Publication des konigl. preussischen 
geodatischen Institutes. Berlin, 1878., 
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Clarke, Col. A. R. Geodesy. Oxford, 1880. 

Fischer, Dr. Philipp. Untersnchuugen iiber die Gestaltder Erde. Darmstadt, 1868. 

Helmert, Dr. F. R. Die mathematischen nnd pbysikalischen Theorieen der lioheren 
Qeodtisie. II Teil. Leipzig, 1884. 

Laplace. M6canique Celeste. Tome 3. Paris, 1878. 

Pratt, John H. A Treatise on Attractions, Laplace's Functions, and tlie Figure of 
the Earth. Fourth edition. London and New York, 1871. 

Stokes, G. G. Mathematical and Physical Papers, Vol. II, especially the paper On 
the Variation of Gravity at the Surface of the Earth. Cambridge University Press, 
1883. 

Thomson and Tait. Treatise on Natural Philosophy. Vol. I, Part II. Cambridge, 
1883. 
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